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rSKRS  MANl'AI. 

I  .  PRINCIPLES  BEHIND  THF.  APPROACH 

Structural  information  is  derived  from  the  distribution  of 
distances  G(r)  (EQ  1)  in  the  sample,  which  is  the  Fourier  sine 
transform  of  the  interatomic  interference  intensity  si(s).  The 
momentum  transfer  s  equals  4rsin0/i  where  X  is  the  wavelength  and 
2e  IS  the  angle  between  the  incident  and  diffracted  beams.  si(s) 
is  not  directly  observable  experimentally.  The  total  diffracted 
intensity  I^  (s)  (EQ  2)  is  composed  of  the  interference  intensity, 

I (s) ,  coherent  and  incoherent  atomic  intensities,  1  (s) ,  plus 

multiple,  I^(s),  and  extraneous,  I^(s),  scattering  terms.  After 
subtracting  the  non-interference  terms,  which  we  collectively  refer 
to  as  the  background  scattering,  R(s)  (EQS  2  and  3),  and  divide  by 
a  sharpening  function,  ueneral ly  i  f  ,  wo  generate  the  interference 
intensity  i(s)  (EQ  4)  corresponding  approximately  to  the  intensity 
from  vibrating  point  atoms. 

To  generate  an  accurate  G(r)  from  a  truncated  set  of  total 
intensity  data,  l^(s),  a  data  reduction  procedure  must : 

1.  Isolate  si(s)  from  I^(s)  by  determing  the  experimental 
background  intensity  H(s). 

2.  Minimize  errors  due  to  termination  and  scaling  of  the 
intensity  data. 

The  radial  distribution  least-squares  program  (RADII..'?)  described 
here  is  characterized  by  a  unique  treatment  of  termination  effects, 
wherein  the  contributions  from  the  first  few  ordered  distances  in 
the  sample  are  subtracted  from  Gi{s)  to  form  si'(s)  (EQ  5),  before 
the  Fourier  transform  which  generates  G'(r)  (EQ  6)  is  performed. 
G'(r)  is  expressed  explicity  (EQ  7)  as  a  function  of  the  total 
intensity,  scattering  factors  and  refinable  parameters  for  the  scale 
and  shape  of  B(s),  the  absolute  intensity  scale  factor,  the  bulk 
density,  and  the  coordination  numbers,  disorder  parameters  and 

interatomic  distances  for  the  first  few  ordered  distances.  Equation 
10  18  minimized  subject  to  the  following  constraints: 

1.  The  inner  region  of  G'(r)  *  zero  where 

distances  are  known  not  to  exist.  By  minimizing  termination  effects 

Mjnu«rnp(  tuhmill^  Aufutt  1979 


which  affect  this  region  as  well  as  the  rest  of  G'(r),  this  region 
can  be  utilized  as  a  sensitive  measure  of  background  error. 

2.  The  outer  region  of  G'(r)  should  show  a  uniform  distribution 
of  distances. 

i.  The  scaling  of  the  intensity  data  should  be  consistent  with 
the  bulk  density  and/or  known  coordination  numbers. 

•J .  H(s)  should  be  a  smoothly  changing  function  whose  shape  is 
compatible  with  theory. 

This  approach  has  several  advantages: 

1.  Terminal  ion  error  is  minimized  because  si'(s)  is  essentially 

terminatovl  at  the  cxperimt>ntal  data  limit  s  .  This  is  the  case 

max 

bec.iuse  the  primary  contributois  beyond  s  are  due  to  the  first 
few  oid»>red  distances  because  of  their  low  disorder  parameters. 

2.  The  structural  parameters  for  the  short  distances  can  be 
accurately  determined  by  minimizing  the  spurious  detail  in  the 
region  of  G'(r)  where  these  distances  contribute.  The  range  over 

ch  G '  ( r )  *■  -I'r*^  IS  recjuired  to  be  zero  is  extended  to  include 

the  region  involving  the  distances  which  are  subtracted. 

3.  Krrors  in  the  structural  parameters  of  the  ordered  distances 
arising  from  the  variability  of  the  scattering  factors  w.th  s,  for 
x-rays  and  electrons,  are  minimized  because  the  short  distance 
contributions,  with  variable  scattering  factors,  are  subtracted 
before  transformation. 

G(r)  IS  generated  by  reintroducing  a  gaussian  description  of 
the  short  distances  without  the  associated  termination  effects. 

Further  details  can  b<*  found  in  the  following  references: 

1.  .J.H.  Konnert  and  ,T .  Karle,  Acta  Cryst  .  (  197  3)  A29 ,  702. 

2.  r.  D'Antonio,  F.  Moore,  J.H.  Konnert  and  J.  Karle,  Trans¬ 
actions  of  the  American  Crystallographic  Association  (1977), 
1_3,  43. 

3.  K.N.  Raymond,  Lecture  Notes,  Least-Squares  Tutorial,  Spring 
Meeting,  American  Crystallographic  Association,  Rerkeley, 

CA,  1974. 
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G(r)  «  47Tr(p(r)  -  p^)  *  —  J  si  (s)  sinsrds 

o 


K(s)  «  I(s)  >  I  (s)  >  I^(s)  ♦  I  (s) 
t  a  rn  e 


(s)  -  I (s)  >  H(s) 


i.(s)  -  (I^(s)  -  H(s))  i.f 


SI 


i(s)  -  si(s)  -  SK^N.  ,,C.,.sin(sr.  ,,)/r.  pxf.(-iJ_s^/2) 


Sd 


AH  AH 


AH  '  AH 


AS' 


sd  refers  to  the  first  few  ordered  short  distances. 


For  electron  diffraction  C,,,  «  f  ,  (s  )  f  „  ( s  )  cos  ( .  ,  (  s )  - 

AH  A  B  A 


,  I  C  !  '  <  ' 

B  '  “  ■ 


for  x-ray  and  neutron  diffraction  •  =  n. 

spax 


G’  (r)  -  4-r  (.  ‘  (r) 


1  )  ■  —  /  SI ' (s) sinsrds 

o 


(f' 


G’(r)  -  SKr.  *  -  (1^3(1  (s)  -  H(s))/  If^ 

°  ■  s 


SKi:N,.^j,C^.^pSinsr^j,exp(-i;,^s‘/2)/r^,Jexp(-.s‘)sinsr-s 


sd 


AH 


AH' 


B(s)  can  be  expressed  two  ways  in  the  pronrarr: 

1.  OVERIJiPPING  EXPONENTIAL  HACKGROUNP  -  single  and  multiple 
ranqe  data 

‘'kn 

B|^(s)  »  Sk5  iRamp(s)  ♦  ^  ( s )  exp  ( A^^^*  ♦  ^km*  '  ' 

r  '  1 

NCV(k)  are  the  number  of  backqround  seaments  in  the  k'^^  backorounc’ 
for  the  overlapping  set  of  data.  For  sincile  ranq«>  data  k  = 


1 


2.  GAUSSIAN  SUMMATION  BACKGROUND  -  single  range  data  only 
NG  2 

B(s)  -  SK5(^'  A  exp(-B.s  )  ♦  CO  >  Ramp(s)) 

1^1  ‘  ^ 

NG  equals  the  number  of  gaussians 
Rampls)  -  KA.s 

Minimize  |  5^  (G'(r)  ‘ 

I  inner 
req  ion 


(10) 
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I  I  .  SCOI’E 


The  Radi a  1  Distribution  Least-Squares  Program  (RADILS)  can  be 
used  to  analyze  amorphous  or  liquid  x-ray,  neutron  or  electron  dif¬ 
fraction  data,  by  providing  the  appropriate  scattering  factors. 

•X-ray  and  neutron  diffraction  scattering  coefficients  are  calculated 
in  RADILS  while  those  for  electron  diffraction  must  be  supplied  from 
a  separate  program.  A  choice  among  several  background  treatments 
and  solutions  of  the  normal  equations  is  provided.  The  structural 
parameters  can  be  constrained  with  linear  constraint  relations  if 
necessary . 

Po lycrys t a  1 1  me  data  can  be  broadened  to  simulate  amorphous 
scattering  for  topology  testing  by  RDK  comparison. 

RADILS  has  been  programmed  to  process  up  to  four  sets  of  over¬ 
lapping  data  covering  different  data  ranqes,  but  this  version  has 
not  been  extensively  tested. 

This  situation  is  common  in  the  sector-microdensitometer  method 
of  electron  diffraction  and  for  x-ray  data  whore  different  slit 
systems  are  used  to  optimize  counting  statistics. 

The  approach  introduces  additional  parameters  into  the  refine¬ 
ment,  i.e.  a  scale  factor  for  each  range  of  data  and  restraints  which 
restrain  the  interference  intensities,  in  the  overlap  regions,  to 
agree  within  specified  limits.  Alternatively,  the  end  points  of  the 
reduced  intensities,  in  the  overlap  region,  can  be  constrained  to 
be  equal.  This  approach  allows  each  data  range  its  own  background 
curve . 
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Ill .  INPUT 

A.  Control  section  (required) 


CONTROL  CARD 

COLS. 

1  (20A4) 

PARAMETKR 

DESCRIPTION 

1-68 

y.TITL(l)  - 

MTITL (17) 

Descriptive  title  card 

69-72 

MTITL (18) 

-  DK43  (Literal) 

Scattering  factors  are  read 

MTITL(18) 

-  blank 

f  rom  file  43 

Scattering  factors  read 

73-76 

MTITL( 19)- 

DK42  (Literal) 

froin  cards 

Experimental  data  are  read 

>!TITL(19) 

-  blank 

from  file  42. 

Experimental  data  are  read 

77-80 

MTITL(20)- 

XRTI  (Literal) 

from  cards 

X-ray  or  neutron  data  are 

MTITL(20)- 

EDIL  (Literal) 

read  in  XRTI  and  processed 

Electron  diffraction  data 

CONTROL  CARD 

1-2 

2  (4012) 

IDX (1) <0 

are  read  in  EDIL  and  processed 

No  DSIS  plot 

3-4 

IDX(2)<0 

No  CSIS  plot 

5-6 

IDX ( 3) <0 

Intensities  will  be  broadened 

7-8 

IDX (4 ) <0 

Print  all  matrices  during 

9-10 

IDX (5) <0 

ref inement 

G'(r)  is  calculated  only  to 

11-12 

IDX (6)<0 

6  A  and  rG(r)  is  not  plotted 

at  all.  Used  in  the  initial 

stages  of  refinement. 

Eliminates  echo  print  of  input 

13-14 

I DX ( 7 ) <  0 

i  ntens 1 1 ies 

No  back  transform  of  AG'(r) 

15-16 

IDX  (8) <0 

is  calculated. 

No  print  of  s,  TI,  SIS,  o(SIS) 

r,  rG(r)  and  o(rG(r))  after 
last  cycle. 
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COLS  . 

rARAMLTKR 

DKSCRIPTION 

17-18 

IDX (9) <0 

Used  for  multiple  data  ranges.  Only 
the  beginning  of  each  overlapping 
reduced  intensity  range  is  constrained 

to  be  the  same. 

IDX  (9)  ■'0 

Both  beginning  and  end  of  each  intensity 

overlap  region  are  constrained  to  be  the 

same.  Not  applicable  for  single  range 

data  NCAS  -  1. 

sO 

1 

o 

IDX(10)<0 

No  constraints  on  the  reduced  intensity 

overlap  regions  at  all.  Not  applicable 

for  single  range  data. 

21-22 

IDX ( 1 1 ) «  0 

SK3(NCAS)  -  the  scale  of  the  last  over¬ 
lapping  reduced  intensity  is  allowed  to 

vary.  Generally  SK3(NCAS)  is  fixed  and 

the  scale  factors  on  all  the  other  ranges 

are  varied.  Not  applicable  for  single 

range  data  NCAS  •  1 . 

23-24 

IDX ( 12) <0 

No  plot  of  TUN, NO  or  n(N,NC) 

25-26 

IDX(13)<0 

Call  LDSIS  and  calculate  distances  for 

a  crystalline  topology. 

27-28 

IDX(14)<0 

Ramp(s)  «  1.0  -  exp(FA.s) 

>0 

Ramp ( s )  -  FA . s 

2<)-30 

IDX (15) • 0 

CSIS  IS  linearly  ext  rafjolated  to  zero 

at  s  »  0 

>0 

TSIS  IS  used  as  inner  region  of  CSIS. 

31-32 

IDX (16) <0 

Klectron  diffraction  data  is  highlighted 

33-34 

IDX ( 17) <0 

AZ,B7.,C7.  and  D7,  are  read  in  for  overlapping 

exponential  bac)cground  in  Cl. 

>0 

Initial  baeV-.ground  points  are  road  in  C2. 

35-36 

IDX ( 18) <0 

Reduced  intensity  is  obtained  by  sub¬ 
tracting  and  dividing  by  the  bac)tground. 

>0 

Bacltground  is  just  subtracted  -  Recommended 

U4 

1 

U4 

GD 

IDX (19) <0 

Punch  CSIS 

39-40 

IDX  (20) <0 

Ramp(s)  -  - 100 . exp (-FA . s) 
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COLS. 
41-42 
4  3-44 


45-46 


47-48 


49-50 


51-52 


CONTROL  CARD 
1-10 
11-20 


21-30 

31-40 


41-50 


51-60 


PARA.MKTKR 
ID.X  (21)  <0 
IDX(22)<0 


IDX(23)<0 


IDX(24)«;0 

>0 


IDX (25) <0 


IDX (26) ^0 


3  (8n0.5) 

SK 

SZ 


DA.Mr 

FA 


SKR 


DAMA 


DESCRIPTION 

Plot  CSIS  and  o (CSIS) 

Eliminates  constraints  on  overlapping 
exponential  background.  Used  only  for 
testing. 

Number  of  cycles  of  smoothing  on  TI  or 
DSIS,  with  square  wave  window  of  2xKRAVG 
{K>  1  n  t  s  . 

Smooth  DSIS  in  three  ranges. 

No  DSIS  smoothing.  TI  is  smooth  if 
KRAVG  ^  0  on  card  4. 

Reduced  intensity  is  obtained  by  sub¬ 
tracting  a  background  composed  of  a  sum 
of  Gaussians.  Alternate  background 
treatment . 

Plot  of  ^*G(r)  contains  individual 

Gaussian  peaks  for  each  of  the  NDIST 
subtracted  and  their  sum.  Gaussian 
decomjxis  i  t  ion  display. 

scale  factor  I  ,  /SK  »  I 
obs 

bulk  density.  This  bulk  density  is  equal 

to  r  ' (  - f  .) ^ /  1 f ^  _  where  r  '  is  the 
o  s  •  0  s»  0  o 

UC  UC 

bulk  density  in  units  of  composition  per 


Experimental  dampinq  factor  exp (-DAMP . s^ ) 
R.mpinq  factor  to  make  intensities  more 
amenable  to  fittinq  with  exponentials. 
Necessary  for  neutron  data. 

If  SKR>0  TI(N,NC)  and  SKRxTUN.NC)  are 
plotted.  Helpful  when  TI(N,NC)  and 
R(N,NC)  are  nearly  equal. 

a  «  -DAMA.r^.  This  is  used  only  in  plotting 
the  final  ^^G(r)  curve  to  reduce 
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61-70 


CONTROL  CARD  4  (615) 


1-5 

NCAS 

6-10 

NSIS2 

11-15 

KRAVG 

16-20 

I  POWER 

the  effect  of  random  errors  at  large  s 
one  feature  of  the  radial  distribution 
curve  at  large  r,  which  arise  primarily 
from  the  low  angle  data.  DAMA  is  read 
as  a  positive  quantity.  If  DAMA  »  0.0 
a  -  DAMP;  if  DA,MA  0  a  *  -DAMP.r^  +  DAMP 
•Maximum  value  of  r  in  the 
■I’^r (r)  -  p^)  plot. 


Number  of  overlapping  data  sets  read  in 
NCAS  “  1  in  present  version  of  program. 
Number  of  intensity  data  included  in 
refinement.  NSIS2  -  (SL£T(NC)  -  s(l)/AS  ♦  1) 
Running  average  parameter.  Window  is 
(2  X  K  RAVC.  *1)  points . 

power  of  r  plotted  in  4’’ r  ^  (o  ( r )  ~  Pq) 

plot . 


B.  Intensity  and  theory  section  (One  of  the  following  two  is 
reifu  1  red )  . 

1.  Kloctron  diffraction  input 

If  MTITL(20)*  P.DIL  Process  electron  diffraction  data. 

If  MTITL(18)»  DF43  Scattering  factors  are  read  from  file  43 

instead  of  cards. 

la.  Klectron  diffraction  scattering  f ictor  coefficients. 

The  electron  diffraction  scattering  coefficients  for 
each  atom  pair  C.„,  FO.  5,  are  generated  in  a  separate 

A  n 

program.  Our  program  F.DCOEFF  is  not  included  in  RADILS. 

The  following  terms  all  apply  to  that  program.  The  user 
m.unt  generate  a  file  for  these  coefficients.  This  approach 
IS  used  because  the  electron  diffraction  scattering  factors 
are  tabulated  at  specific  s  values  and  all  of  the  calculations 
to  generate  C,„  need  not  be  repeated  each  time  RADILS  is  run. 

Ad 
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(2I5,4X,Z4,4X,2A4) 


1-5 

.VCOEKF 

number  of  scattering  factor  coefficients 

6-10 

NSIS2T 

numb*-r  of  scattering  factor  coefficient 
fJCint  s  . 

ii-i ; 

SKIP 

15-18 

S’Ame  1 

Optional  name  of  job  which  created 
scattering  factor  coefficient  file. 

l<»-22 

SKIP 

23-26 

NA.M1.  2 

Optional  date  of  job  which  created 

27-  30 

NA.MK  3 

scattering  factor  file. 

(20A4J  (I  . 

1 . VCOEKF  Road 

the  next  series  of  cards  including  the  END 

card  NOOEKK  times. 

1-4 

1 

fir.st  part  of  atom  pair  coefficient  name 
eg  Si 

5-fl 

NA-Mt:  2 

dash  ( - ) 

0-12 

'.  A.y»:  3 

Second  part  of  atom  pair  coefficient  name 
e.g.  0.  Complete  name  is  Si  -  0 

1 8K10. 5) 

I-IO 

n'xi 

not  used 

11-20 

TYPE 

code  for  type  of  sharpening  function  used, 

O 

1 

wave;.' 

wavenumber  2V  or  4"/\ 

31-4 : 

.SM'; 

minimum  s  v.ilue 

41-SO 

SM.X 

maximum  s  value 

51-60 

facto P 

not  used 

2r 10. 5, 153 

1-10 

s(l) 

first  s  value 

1 1-20 

DS 

s  interval 

21-25 

NSIS2T 

number  of  points 

{5E14. 7) 

ERARf  7.  1) , J  -  1 

,NSIS2T  -  scattering  factor  coefficients 
for  the  atom  pair 

(20A4  ) 

1-4 

LEND 

END  card 
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Read  in  cl«*ctron  diffraction  sharpen  inij  factor 
(2  0A-S) 

1-80  NTITL  "SHAKr^NING  FACTOR"  title 

(8F10. 5) 

1-80  DUMl ,TVrF,KAVKSO,SMN ,SHX , FACTOR  Same  as  for  scatterinq 

factor  coefficient 

(2F10. S, IS) 

1-2S  3(1) ,DS,NSIS2T 

»  (Si:l4.7) 

(FC(I),I-1,  NSI.S2T)  Sharpening  function 

(20A-1) 

1-.J  I.KND  end  card. 

Read  in  electron  diffraction  total  atomic  intensity 
(20A.1) 

1-80  STITI,  "TOTAL  ATOMIC  INTF.NSITY"  title 

(8F10. S) 

1-80  Dl'Ml  .TYFF  .WAVF.NO,  SMN  ,SMX,  FACTOR  Same  as  for 

(2F10.S,IS)  scattering  factor 

1-25  s(l),  DS,  NSIS2T  coefficients 

(5KU  .  7) 

(T1 A  (N)  ,  N' 1 ,  S'S  IS2T)  Total  shariK*ned  atomic  intensity 

(20A4  ^ 

1-4  LFS'O  card. 

lb.  Klectron  di  f  f  r.act  ion  leveled  intensity  if  MTITL(19)  «  DF42 
the  intensity  is  read  from  file  42  instead  of  cards.  This 
intensity  is  the  output  from  our  leveling  and  averaging 
program  GFTIL.  RADII.S  ex{>ects  intensity  data  that  has  been 
correctevl  for  systematic  errors  and  sharpened.  If  NCAS  >1, 
read  this  set  of  cards,  including  FNO  card,  NCAS  times.  The 
pres«*nt  version  has  only  been  tcstetl  for  NCAS  •  1  single 
range  data. 

120A4) 


1-8 


Title  ca  rd 


1 


(8K10. 5) 

1-10 

CH 

Camera  heiqht  in  nm 

11-20 

I'TOS 

plate  to  rotating  sector  separation 

21-30 

WAVENl) 

wavenumber  4w/X 

31-40 

SKRST(NC) 

first  data  print  for  the  NC^^  data  range 

41-50 

SL5T(NC) 

last  data  point  in  the  NC^^  data  range 

51-60 

FACTOR 

not  used 

61-70 

SF(NC) 

(SK/SK 3 (NO )  necessary  to  bring  the 

theoretical  atomic  intensity  to  the  level 
of  the  NC*^*'  data  set  for  subtraction  in 
KDIL.  If  SF(NC)  •  0.0,  no  atomic  theory 

IS  subtracted. 

71-00 

PCTSG 

When  standard  deviations  from  counting 

statistics  are  not  available  as  in 

electron  diffraction  photographic  data 

the  error  in  the  intensity  is  expressed 

as  percent  error  envelope.  The  error 

IS  expressed  as  a  p>erccnt  of  the  sharpening 

function.  PCTSG  “  0.005  corresjxinds 

to  0.5»  of  the  sharpening  function. 

(2F10. 5, 15) 

1-10 

s(l) 

first  value  of  s 

11-20 

DS 

s  interval 

21-25 

NSISl 

number  of  points 

(5K14 . 7) 

(TI (N,NC) , I-l . 

Total  intensity  data  for  the  NC 

NSISl ) 

data  range 

(20A4) 

1-4 

I,FND 

F.ND  card 

12 


L 


1 


2.  X-ray  and  neutron  diffraction  input 

If  M"ITL(2Q)“  XRTl  Process  x-ray  or  neutron  diffraction  data. 
If  MTITL(19)-  DK42  Read  data  from  file  42. 

2a.  Kxperimental  intensities. 

If  multiple  overlappinu  sets  of  data  are  beinq  processed 
(NCAS  •»  1)  read  the  followinq  set  of  cards  NCAS  times. 

Present  version  tested  for  NCAS  «  1  only. 

(1615) 


1-5 

NSISl 

number  of  total  intensity  data  points 

(9P10. L. 

1-10 

s(l) 

f irst  value  of  s 

11-20 

DS 

s  interval 

21-30 

SFRST (NC) 

first  data  p>oint  in  the  set  of  data 

31-40 

SLST(NC) 

last  data  p>oint  in  the  set  of  data 

41-50 

SCTI 

scale  factor  modifyino  TI(S,NC)  and 

SG(N,NC)  before  refinement.  Usually 

set  to  1.0 

(20A4) 

1-00 

NFMT 

variable  input  format  for  total  intensities 

ea  ( (4 (2F10. 5) ) ) 

(20A4) 

1-80 

S7I7L 

Title  card 

(SFMT) 

2b. 

(TI (N,NC) ,SG (N,NC) ,N-1 , NSISl )  7otal  intensity  and  standard 

deviation  of  the  NCth  data 
set 

Scattering  factors 

If  M7ITL(18)  «  DF4.'>  -ead  theory  from  file  43  instead  of 

cards. 

(9F8. 4 ,F5 

. 2 )  M- 1 ,  NDA 

Bead  one  card  for  each  of  the  NDA 

distinct  atom  types. 

1-8 

AA(I,M) 

The  AA,  BB  and  CC ' s  are 

9-16 

BB(1 ,M) 

the  coefficients 

17-.''4 

.\A  t  2 

M  1 

tor  analytic  approximations 

2 ‘.-3  2 

Bl' t  .■ 

V) 

of  Haitreo-Fock 

X-ray 

3  1  -  i  I 

AA  3 

M) 

scatter inq  factors.  Values  are 

41-46 

BH(  ’ 

'*) 

tabulated  by  B. 

Mann,  (Acta 

49-St 

AA  4 

M) 

Cryst.  (1968), 

A24 ,  32)) 

"  '-^4 

'  4 

M  ' 

r,  J 

00  ,  V 

•  1  _  :  ■ 

Ml  v.v; 

1 

number  of  atoms 

in  unit  of  composition 

'8-fln 

-ar.l 

label  in  1 

information,  not  read 

t  e  1 

-  K'.  1  neutron  dif 

fraction  data  AA, 

PB  and  CC  are  zero 

and  OU(V)  is  the  coherent  neutron  scat  ter inq  amplitude 
tabulated  by  0.  K.  bacon  (Acta  Cryst .  (1972)  A^fl.  357). 


The  pair  typo  indicator  on  the  short  distance  parameter 
■card  Nr(b)  IS  dotormino^l  by  the  order  of  the  scat  ter  inq 
:act',>r  cards.  Kor  example  it  the  scatterinq  factor  for 
b  1  IS  first,  followed  by  the  O  card  for  SiOj 

sr(I.}  DISTASCr 

1  .91  -  Si 

2  .Si  -  O 

j  .  -  ,1 


Ini’ial  I' I  li  around  Input  (1  of  the  followinu  three  sections 
■  'i  t  e  :  11  r<‘  t'  . 

1.  If  IIiXil"  n  re  1  1  ’’v.'r  1  app i n<j  Exponential  rarameters. 

T);e  I  acV.  ::  'und  is  calculated  directly  from  the  expression 
liven  in  Fc. .  8  Tresent  version  is  only  valid  for  NCAS  ■  1. 


IlO.iFlO.i! 

•iC'  1  ,nca;: 
I - 1 0  NOV ( NO ) 

11 -20  .SK5(NC) 

21-30  SK3(NC) 


If  NCAS  «  1  read  NCA.s  cards 
number  of  bacicqround  seqmcnts  in  the 
ranqe  of  overlappinq  data 
(Overall  bac)tqround  scale  factor  for 
the  Nc''^  ranqe  of  overlappinq  data 
scale  factor  scalinq  experiment  to  theory 
for  the  ranqe  of  data.  Usually  set 

to  1.0  and  refine  SK. 
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il-40  SFRSTINC)  First  experimental  {>oint  for  the 

data  set . 

41-SO  SLST(NC)  last  experimental  point  for  the 

data  set  included  in  the  refinement, 
NSIS2  on  card  4  must  be  consistent 
with  SI^T(NC) 

NSIS2  -  (SI.ST(NC)  -  s(l)/DS  +  1) 

Nt*V(N'"  are  the  number  of  background  se<7ments  in  the  back¬ 

ground.  If  NCAS  sets  of  overlapping  data  are  being  refined 
there  will  be  NVAS  sets  of  NCV(NC)  background  segment  cards. 


This 

Version  IS  tested  for 

NCAS  •  1  only. 

4K14 . 7 1 

1-14 

NC*  1  ,  SCA.^' 

AZ(1,NC) 

OV'er  lapping  exponential 

15-29 

HZ  a  .NO 

background  segments 

1 

W 

CZ l 1 ,NC) 

for  the  NC^^*  background. 

1  1-46 

DZ  (1  ,NO 

If  multiple  overlapping  data  ure 

1-14 

AZ  tNCV(NC)  , 

NO 

refined  read  this  set  for 

each  range*.  Present  version 

NC  «  1  only. 

15-29 

HZ (NCV(NC) , 

NO 

29-  12 

CZ  (N('V<NC)  , 

NO 

3  1-46 

0Z(NCV..'4C)  , 

NO 

;16F5.2' 

1-5 

NO  1  ,  NCAS 

.SNN  1  1  ,NC) 

These  are  the  points  at  which  the 

6-10 

SSN (2 ,NC) 

overlapping  exponential  background 

11-15 

SSNI l.NC) 

segments  are  constrained  to  be 

• 

SSN (NCV(NC) 

,NC) 

equal  with  I-igranqian 

multipliers.  The  points  are 
generally  set  equal  to  the 

midpoints  of  the  background 
segments  of  length  As.  Road  one  card 
for  each  overlapping  set  of  data. 
Present  version  NCAS  ■  1  only. 
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Kxainple  -  for  3  ovorlappinq  segments  of  length  As  ■  4 


0  2  4  6  8 

S 

SNN  ( 1  ,  1 )  •  2.0  seqments  1  and  2  are  const  r.»  mod  equal. 

SNN(2,1)  •  4.0  seejments  1  .ind  3  are  constrained  equal. 

SNN(3,1)  "  fi .  0  seuTnents  2  and  3  are  constrained  equal. 

2.  ID.X(17)  ;  0.  H<>ad  oveilappini  i‘xt>on«-nt  la  1  trial  {xiints 

This  IS  the  recommended  background  start  inn  point. 

Present  version  only  tested  for  NCAS  »  1. 

{I10,7K10.4)  NC-l.NCA.S  If  NCAS  •  1  read  NCAS  cards 

1-10  NCV(NC)  numlier  of  background  seuments  in 

the  ot'erlappina  set  of  data. 

11-20  SKS(NC)  overall  background  scale  factor 

for  the  Nc'’^  overlapping  set  of  data. 

21-10  .SKI  (NO  scale  factor  placing  the  Nc''^ 

reduced  intensity  on  an  absolute 
scal<’.  Usually  set  to  1.0  and  SK 
ref ined . 

31-40  SKRST(NC)  first  ex|>er  imenta  1  point  in  the 

NC*^^  data  sot. 

41-50  S1.ST(NC)  last  experimental  {.xiint  in  the  Nc'^^ 

data  set.  NSIS2  on  card  4  equals 
(SLST(NC)  -  s(l) )/DS  ♦  1)  . 

(8F10.4)  NC-1,NCA.S  If  NCAS  >  1  read  NCAS  cards 

1-10  SSd.NC)  SS(N,NO  and  FF(N,NC)  are  the 

11-20  FFd.NCI  s  value  and  trial  background  value, 

th 

21-30  SS(2,NC)  respectively,  for  the  NC  overlapping 

31-40  FF(2,NC)  set  of  data.  The  background  parameters 

.  .  are  determined  by  a  three  point  fit 


If) 


.  of  successive  trial  jxiints 

eg  123;  234;  345  et. 

SS(NCV(NC)  ♦^2,NC) 

FK(NCV(NC)  4-2,  NO 


Note  1.  The  selection  of  the  SS(N,NC)  are  based  on  the  segment 

1st  1st 

length  As  which  should  be  >  2i/r  .  r  is  the 

first  bonded  distance  expected  in  sample.  This  is 

important  early  in  the  analysis  to  prevent  background 

errors  from  interfering  with  the  structural  parameters. 

As  need  not  be  constant  and  in  fact  one  or  two  segments 
1st 

<  2',  r  in  the  inner  region  (less  than  4s)  can  be 

helpful  in  providing  the  curvature  needed  to  fit  the 

experimental  background.  As(N)  *  SS  (N  ♦  2, NO  -  SS(N,NO 

Note  2.  If  FF(l,NO  is  zero  then  the  initial  background  is 

determined  automatically  from  a  fit  to  the  theoretical 
2 

If  which  IS  ramped  and  scaled  to  generate  an  si(s) 
which  has  e<jual  positive  and  negative  area.  This  is 
accomplished  by  minimizing  the  following  relation 
N.SISl 

E  s(S)/FC  (S)  (TI (N,NO - (C, *F0  (N)  ♦  C,  ♦  C  ‘s(N))l^ 

N-1 

this  starting  point  may  be  helpful  for  x-ray  and  neutron 
data 


(16F5.2)  NC-1,NCAS  If  NCAS  '  1  read  NCA.S  cards 

1-5  SNN(l,NO  These  are  the  points  where 

6-10  SNN(2,NO  the  overlapping  exponential 

11-15  .  segments  are  constrained 

16-20  .  to  be  equal .  Example 

.  .  in  Cl  and  below. 

SNN  (NCV(NO  ,NC) 
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1 

Tl 


Ji'll  wu'.'i  i-x-jr-r  li-  P-..iy  hi-lp  in  d<‘tprmininvi  .‘.s, 
» *.  ,KK  and  rsNiN.Nt'l.  if  the  sample 


N 

KK  ( ,  Nf  ) 

.SNt;  (N 

1 

0.05 

1.0 

2.0 

■» 

2.0 

0.  H8 

4 . 0 

1 

4.0 

0. 58 

6.0 

•« 

t.  .0 

0 

0 

5 

CD 

0 

0.15 

I.  Ili.\  a.  hf.i.i  ''■.,i;;ss  i  wi  >  ipjri.iti  )n  I'.irmot  ox  s 

TheS'-  initi  il  t. It. «?•■•••  ors  .III-  lot  or  m  i  n«'ij  from  a  separate 
:i  ).:r,ir.  ;a:  SHh  ty  :;;*i:iJ  the  It  .  This  typo  ol  background 

txoatmor.t  .>(5  lios  .'nly  *0  smile  ran  lo  data  NCAS  ^  1. 

r.ist  p.o  sv  that  it/vVi  Is)  KA.s  The  backuiound 

.s  loioiminod  tiy  ’ho  oxj  t  ess  ion  fiivon  in  K'j.  '♦  • 

I  ,  ^X,  il  l''  .  ' 

*>  number  of  la  us  si. ins 


.-;o  ssif 

ll-.’C  KA 

:  1  -  I 

1 1  -  4  . .  b  h  ')  I  1  ; 


r.imp  par.iPW't  or 
.-ons*  ant 

uvorall  back around  scale  factor 


8FI.T.  5 

l-l.l 

AX  ;  1  .  1 

1  1-20 

BT  '  1  ,  1 

21-  1. 

AZ  (2,1 

11-40 

BZ  (  2  ,  1 

, 

AXi!»,l  and  B7.  \N,l!  are  the 
la.issian  amplitude  and 
t  roath  p.irameters,  respect - 
i  ve 1 y . 


A7.  .  N<;,  1  ' 

BZ  isr.,  1  i 


D.  IDX(3)  <  0  BROAD  input  (optional) 

Used  to  broaden  a  po lycrys t a  1 1 1 ne  pattern 


(8F10. 5) 

1-10 

I’SZ 

radius  of  the  ordered  reqion. 

11-20 

SLOl’E 

slope  of  a  linear  bac)turour'.d  which  is  subtracter: 
from  total  int«*nsity  liefort'  liroaden  i  nu  . 

21-30 

XI  AM 

wavelenqth 

31-40 

BINT 

intercept  of  linear  bac)cciround 

K.  IDX(13)  <  0  LDSIS  input  (optional) 

Coneratfs  distances,  coordination  numbers  and  thermal  parameters 
for  a  crystalline  tofxiloqy. 


015) 


1-5 

NATTY 

♦'-10 

NATAS 

11-15 

NSYV 

16-20 

NXl 

21-25 

NX2 

26-  30 

NYl 

31-  35 

SV2 

36-  40 

NZl 

41-45 

.NZ2 

(8F10. 5) 

1-10 

DISMN 

11-20 

DISMX 

21-30 

AA 

31-40 

BB 

41-50 

CC 

51-60 

CA 

61-70 

CB 

71-80 

CG 

namb«T  of  distinct  atom  typ'es 

numl>er  of  atoms  in  the  asymmetric  unit 

number  of  s^-mmetry  t  ransfornat  ion  cards 

defines  the  number  of  unit  cells  to 
search  in  x  direction  foi  contri- 
butinu  atoms 

defines  the  number  of  unit  cells  to 
search  in  the  y  dirt'ction  for 
contributinu  atoms 

defines  the  number  of  unit  cells  to 
search  in  the  z  direction  for 
contributinu  atom.s 

Minimum  distance  of  calculation 
Maximum  distance  of  calculation 
these  are  the  unit  cells  dimensions  and 
direction  cosines. 


10-DISV.X,  AA 
10*DISMX  'AA 

lO-DISMX/BB 

lO+DISMX/BH 

10-DISMX/CC 

10+DISMX/CC 
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(  3F10 . 5 , 1 1 0 , F 10 . 5 )  Read  one  card  for  each  atom  in  asymmetric 

unit  N-1,NATAS. 


1-10 

.XX  (N) 

Fractional  coordinates  of  the  atom 

11-20 

YY(N) 

in  the  asymmetric  unit. 

21-30 

ZZ  (N) 

31-40 

NNTYPF.  (N) 

atom  type 

41-50 

WT(N) 

Weight  used  to  generate  proper  coordination 

nv.imber.  The  weight  can  be  set  to  be  the 
reciprocal  of  the  square  root  of  the  number 
of  atoms  of  the  l^*"  type,  NTyPE(l),  in  the 


asymmetric  unit. 


(12F5. 2) 

Read  one  card 

for  each  symmetry  operation. 

NSYM  cards 

1-5 

DX(N) 

DX(N) 

6-10 

FXl (N) 

DY(N) 

IS  the 

¥  K 

'  translational  vector  for  the 

11-15 

FX2 (N) 

DZ  (N) 

N  symmetry 

ope»  ation 

16-20 

FX3  (.V) 

21-25 

DY(N) 

FXl (N) 

1  FX2(N) 

FX3 (N) 

is  the  rotational  symmetry 

26-30 

FYl (N) 

FYl (SI 

1  FY2(N) 

FYl  (N) 

transformation  for  the 

31-35 

FY2(N) 

FZl (N) 

1  FZ2(N) 

FZ3 (N) 

symmetry  operation.  See 

36-40 

FY3(S) 

symmetry  transformation  in 

41-45 

DZ  (N) 

Glossary . 

46-50 

FZl (N) 

51-55 

FZ2  (.S) 

56-60 

FZ3 (N) 

(8F10. 5) 


1-10 

RDD 

value  of  the 

first  bonded  distance 

11-20 

WOO 

amplitude  for 

first  bonded  distance 

21-30 

RDD2 

distance  and 

amplitude  of  any  other 

31-40 

WDD2 

distance  for 

which  a  good  approximation 

of  the 

amplitude  is 

known . 

41-50 

A.MPMAX 

maximum  val'e 

of  the  amplitude  at  large 

values 

of  r . 
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F.  Short  distance  parameter  imput  (required) 


(1615) 

5  NDIST 

(F10.6, I10,2F10.6) 

I- 10  A(L) 

II- 20  NP(L) 

21-30  X12VB 

31-40  RZ(L) 


number  of  short  distances 

Read  NDIST  cards.  One  for  each  short  distance 
L-1 , NDIST 

Coordination  number  for  the  distance 

Pair  type  of  the  distance.  Determined 

by  the  position  of  the  scattering  factor  cards 
in  Bla  and  B2b . 

B(I.)  -  -XIJVB**2/2.0  .  XI2VB  is  the  root-mean 
square  amplitude  of  vibration  for  the  current 
distance . 

The  interatomic  distance. 


Note:  If  IDX(7)  <  0  G(r)  will  be  bac)c  transformed  out  to  RZ(NDIST) 

to  get  a  better  approximation  of  the  overlapping  exponential 
background  p>oints.  These  pxDints  are  printed  after  the  total 
intensity-background  plot  and  punched  on  unit  41.  They  can 
be  used  as  starting  points  when  resubmitting  the  job. 


G.  Refinement  section  (required) 

The  KI (N)  vector  specifies  which  parameters  are  refined  and  the 
type  of  background  treatment.  There  are  a  total  of  NPAR  flags. 
NPAR  •  3*NDIST  *  4.  The  first  two  flags  refer  to  scale  factor 
and  bulk  density,  the  next  3*NDIST  refer  to  the  structural 
parameters  of  the  NDIST  short  distances  and  the  last  two  refer  to 
the  type  of  background  treatment  and  the  scaling  of  the  observed 
data.  Generally  the  last  flag  is  set  to  1  and  SK  the  scaling  of 
theory  to  experiment  is  varied.  The  background  flag  is 
NPAR.M  »  NPAR  -1  and  the  structural  parameter  flags  end  at 
NNP  NPAR  -  2. 

(512,3512, (14012) ) 

1-2  IDl  IDl  "  1  eigenvector-eigenvalue  procedure 

is  used  to  obtain  parameter  shifts. 
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IDl  “  0  matrix  inversion  is  used 

3--J  NCOS’  number  of  imposed  constraints  on  structural 

parameters.  Background  constraints  are 
applied  automat ica 1 1 y  in  proqram  using 
Laqranqian  multipliers. 

5-6  ID3  not  used. 

7-6  NADD  NADD  «  1  Gaussian  description  of  short 

distance  peaks  are  reintroduced  into  the 

•It  (.'  (r )  -  ^  plot . 

SADD  *  0  short  distances  are  not  reintroduced. 

9-10  SSUB  NSUB  »  1  Plots  of  TI(S,?JC),  n(N,SC),  SIS(N,NC) 

CSIS(N),  o(CSIS(N)),  G(r),  (r)  - 

are  displayed  before  current  refinement  cycle. 

11-12  KI(1)  KI(l)  “2  SK  scale  factor  is  refined. 

KI(1)  =  1  SK  scale  factor  scalinq  theory  to 
to  experiment  is  not  refined. 


13-14 

KI  (2) 

KI  (2) 

■  3 

S2  bulk  density  is 

ref i nod . 

KI  (2) 

-  1 

SZ  bulk  density  is 

not 

ref ined . 

15-16 

KI  (3) 

KI  (3) 

-  4 

Coordination  number 

of 

first  short 

distance  is  refinf^d 

. 

KI  (3) 

-  24 

Coordination  number 

of 

f  1  rst  short 

17-18  KK-J) 


distance  is  constrained.  Q  imatrix 
will  specify  which  other  coordination 
numbers  are  involved  in  the  constraint 
relation. 

KI(3)  =  1  Coordination  number  is  not  refined. 

KI(4)  -  5  The  first  short  distance  is  refined. 

KI(4)  «  25  The  first  short  distance  is  constrained. 

0  matrix  will  specify  which  other  dis¬ 
tances  are  involved  in  the  constraint 
relation. 

KI(4)  »  1  The  first  short  distance  is  not  refined. 
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19-20  KI(5) 


21-22  KI(6) 

23-24  KI(7) 

25-26  KKfl) 


KI(5)  “  6  Thermal  parameter  for  the  first  short 
distance  is  refined. 

KI(5)  «  26  Thermal  parameter  for  the  first  short 
distance  is  constrained.  The  Q  matrix 
will  specify  which  other  thermal 
parameters  are  involved  in  constraint 
re  la  t ion . 

Coordination  number  flag  for  second  snort 

distance 

Second  short  distance  flag. 

Thermal  parameter  flag  of  second  short  distance. 

The  same  numbering  scheme  as  in  KI(3),  Kl(4), 

and  K I (S)  is  used . 


KI  (S'N'P)  Thermal  parameter  of  the  NDIST  distance 


KI(NPARM)  KI(NrAR.v.) 

-  7 

A  and  h  of  each  overlapping 

exuonential  background  segment  is 

ref ined. 

KI (NPARM) 

■  8 

A  and  C  of  each  overlapping 

exponential  background  segment 

is  refined. 

KI (SPARM) 

-  9 

SK5  overall  scale  factor  of  over¬ 
lapping  exponential  background  is 

ref Ined . 

KI (NPARM) 

-  10 

A,B  and  C  of  each  overlapping  ex¬ 
ponential  background  segment  is 

ref ined. 
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KI (NPARM) 


12  A,B  and  D  of  each  overlapping 
exponential  background  seament 


IS  refined. 

KI (NPARM)  ■  13  SK5,  CO,  FA  of  Gaussian  summation 
background  is  refined. 

KI (NPARM)  •  M  A  and  B  of  each  qaussian  is 
ref ined . 

KI (NPAR)  KI (NPAR)  •  11  SK3  scaling  of  experiment  to 

theory  is  refined. 

Usually  set  to  1.0  and  not  refined 

KI(NPAR)«  1  .SK3  not  refined. 

It  NCON  •  0  Re.id  Q(N,M)  matrix  which  specifies  which  of  the 

original  nfARR  parameters  read  on  the  KI  card  are  to  be  varied 
independently  (uppvcr  section  SIV  x  SPARR)  and  how  the  dependent 
parameters  are  related  (lower  section  NCON  x  NPARR)  .  The  uppicr 
section  IS  automatically  determined  in  the  program  from  the  KI 
card  i.e.  all  flags  less  than  ^  20.  The  lower  section  is  read 
from  NCON  data  cards.  These  cards  represent  the  last  NCON  rows 
of  0-  Each  row  is  compxDsed  of  the  coefficients  of  all  parameters 
involved  in  the  linear  constraint  relation.  Read  one  card  for 
each  constrai-'.t  relation.  NCON  cards. 

(10F8.3)  L  -  NPARR  -  NCON  ♦  1 

1-9  0(L,1)  If  these  parameters  are  involved  in  the 

9-16  0(b,2)  constraint  relation,  these  are  the 


17-24 


0(L,  3) 


coef f icients. 


KXAMl’LK  OK  KI  CARD 


Type  of  Refinement 

N 

KI  (N) 

Parameter  Refined 

1 

2 

SK  (scaling  of  theory 

to  experiment) 

2 

3 

SZ 

3 

4 

SI 

4 

5 

rl 

Structural  Parameters 

5 

(. 

11 

Kq,  9 

6 

24*  * 

N2 

Parameters  can  be 

7 

5 

rl 

constrained. 

8 

6 

12 

9 

1 

not  refined  N(NDIST) 

10 

5 

r (NDIST) 

N’ST 

6 

1 (NDIST) 

Overlapping 

NPARV* 

7 

A,B  (Suggested  paramc 

Exponent  1  a  1 

8 

A,C 

Background  Eq.  10 

9 

SK5 

10 

A,B,C 

12 

A,B,D 

Gaussian  summation 

1  3 

SK.CO.FA 

Background  Eq .  1 1 

14 

A,  B 

NPAR 

11 

SK3  (scaling  of  exper 

— 

to  theory) 

•  CHOOSE  ONLY  ONE.  If  KUNPAR,”)  - 

1  skip 

background  refinement 

**  Tho  second  coordination  number  is  constrained.  0  matrix  specifies 

which  other  coordination  numbers  are  involved  in  constraint  relation. 
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.2.2 (F5.2, I5,2rS.2) ,3F5.2) 

1-5  SFT  Applied  shifts  are  SFT  times  calculated  shifts 

6-10  FOBS(l)  first  r  value  in  the  first  observational  range 

of  the  least  squares.  This  is  the  inner  region 
of  G'  (r)  ♦  (•>’'rr(r). 

11-15  NOBS(l)  Number  of  observations  in  the  inner  region  of 

•In  1(1  (r) 

16-20  DOBS(l)  increment  in  r  in  the  inner  region  of  4'-r,  (r> 

21-25  WOBS  ( I )  weight  of  the  inner  region  of  -I'r.Mr)  in  the 

ref inement . 

26-30  FOBS (2)  First  r  value  in  the  second  observational  range 

of  the  least  squares.  This  is  the  outer  region 
of  G' (r ) . 

31-35  NOBS (2)  Number  of  observations  in  outer  region  of 

G'(r).  Set  to  1  if  nist  refining  inner  region. 

36-40  DOBS (2)  increment  in  r  in  the  outer  region  of  G'(r). 

41-45  WOBS(2)  weight  of  the  outer  region  of  G'(r)  in  the 

ref inenent . 

46-50  WOBS(3)  weight  of  the  third  observat lona 1  range. 

51-55  WOBS (4)  weight  of  the  fourth  observational  range. 

56-60  WOBS(5)  wt'ight  of  the  fifth  observational  range. 

These  last  three  observational  ranges  are  used 
when  overlapping  sets  of  data  arc  being  refined. 
They  ari>  the  three  possible  overlap  regions  in 
the  reducetJ  inten.sity  from  a  maximum  of  four 
overlapping  sets  of  total  intensity  data. 

The  retiuced  intensity  CSIS(N)  must  be  continuous 
and  therefore  the  twc  intensity  contributors 
in  an  overlap  region  must  be  restricted 
This  can  be  done  with  constraints  on  the  end 
points  or  by  restraining  their  difference  to  be 
below  a  certain  limit.  These  weights  are  used 
when  the  overlap  regions  are  restrained.  The 
constraints  are  selected  with  the  IDX  control 
card.  If  constraining  these  weights  are  0. 
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I 


If  IDl  >  0 
{3F10. 5, 110) 

I-IO  EVMIN 

11-20  RSMIN 

21-30  SHIFT 

31-40  MAXPAR 


minimum  eiqenvalue  accepted  in  filtering 
process  in  IGN.  If  0  all  eigenvalues  are 
accepted . 

minimum  residual  reduction  required  in 
filtering  process  in  IGN  for  a  particular 
eiqenshi f  t . 

maximum  eioenshift  (linear  parameter  shift) 
allowed  in  filtering  process  in  IGN. 
maximum  numijer  of  eigenvalues  (linear 
parameters)  all'-wed  in  filtering  process 
in  IGN.  Program  selects  the  MAXPAR  eiqen- 
shifts  which  produce  the  greatest  residual 


rei'uct  ion . 


IV.  OPERATIONAL  REMARKS 
A .  Constraints 

The  interatomic  distances  ;rZ(L)),  the  coordination  numbers 
(A(L)),  and  the  disorder  parameter  (B(L))  can  be  constrained  by 
NCON  linear  constraint  relations,  0,  which  are  read  in  RADILS. 

If  NCON  »  0  no  constraints  are  used.  me  imposition  of  constraints 
lessens  the  number  of  variables  (NPARR)  and  required  suitable 
combinations  of  the  oriqinal  derivatives  i.n  constructing  a  new 
derivative  matrix  for  the  smaller  set  of  independent  variables, 

NIV.  The  new  derivatives,  with  respect  to  the  independent  para¬ 
meters,  V,  can  be  obtained  from  the  old  derivatives,  with  respect 
to  the  original  dependent  parameters,  x,  by 


( 


dF 


where  F  is  the  function  being  minimized,  T  is  the  transpose  operation 
and  J  •  0  .  The  original  parameter  chamies  are  obtained  in  RESET 

from  the  independent  set  by  ux  ■  J  .W  . 


The  0  matrix  consists  of  an  upper  sectuon  U  (NIV  x  NPARR)  and 
a  lower  section  L  (NCON  x  NPARR).  The  upper  section  determines 
which  of  the  original  parameters  will  bo  treated  as  linearly  inde¬ 
pendent  variables  and  the  lower  section  su-mmarizes  all  the  constraints 
imposed  on  the  problem,  by  specifying  the  relationships  between  the 
dependent  and  independent  variables. 


y 

dv 

dx  •  j 

L 

0 

0  dx  •  dv 


The  upper  section  is  automatically  calculated  from  the  KI  card. 
When  a  structural  parameter  flag  4,  5  or  6  is  entered  as  24,  25  or 
26  that  parameter  is  considered  to  be  constrained.  U  is  found  by 
placing  ones  and  zeros  in  appropriate  locations. 
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Thf  low»>r  soi't  ion  L  is  read  on  the  cards  followinq  the  KI  cara 
\JN  •  0.  The  coefficients  of  the  linear  constraint  relations 
laced  hy  the  user  in  the  appropriate  column.  The  record  length 
•ach  constraint  relation  is  NPARR  x  8  fields.  If  NPARR  is  >  10 

• han  one  card  will  be  required. 

Txam; le:  Silica  ulass  -  We  would  like  to  vary  the  scale  factor 

Oik  l.-nsity  SZ,  the  coordination  numbers  of  the  Si-O,  O.  .  .0  and 
.  and  the  .  verall  background  scale  factor  SK5. 

In  ad  Ilf  ion  we  would  like  to  constrain  the  coordination  numbers 
1  1 '  ws 

.2  .  ..1 


l.S  X 


;  *  — w 

'  S 1 ■  O.  .  .0 


N 


Si  -O 


•’Si-'. 


.3 

■’Si .  .  .Si 


O.  .  .0 


S  =  N 


Si . . .Si 


■d . at  icn  cl  Normal  Fquation 

Th«'  user  has  a  choice  in  obtaining  the  shift  vector.  Either 
'c.TFT  M  IdN  ci:-.  be  used  dependinu  on  the  complexity  of  the  problem. 
i  Il.'.  -l  I  'c;  IS  selected  and  the  shifts  are  determined  by  the  cigen- 
1 1  1  lenvect  .ir  technique.  Several  limiting  parameters,  read  in 

,  can  help  wi^h  unstable  matrices  by  filtering  out  eigenvalues 
nicb  ar--  either  too  small,  cause  excessive  shifts  in  the  linear 
■mbina'i  'r.  parameters  or  insufficient  reduction  in  the  linear 
;;roximitinn  residual.  Since  the  eigenshifts  are  uncorrelated, 

■•ting  any  of  them  to  zero,  if  it  does  not  meet  a  particular  criterion 
=  s  net  affect  the  required  values  of  the  remaining  ones.  Thus 
ilterinu  provi  les  a  means  of  ensuring  that  large  eigenshifts  will 
■  t  be  illowe-l  unless  the  resultant  residual  reduction  is  judged 
.tficient.  It  IS  even  possible  to  refine  correlated  parameters 
imvi  1 1  anecuR  1  y ,  using  1(3N  i.e.  SK  ,  SZ,  N,  R1  ,  SK5. 


Invert  uses  the  ASC  system  matrix  inversion  DSIMI'C,  which  is 
a  double  precision  routine  operatinq  on  a  linear  packed  array 
APACK.  SIMEC  (a  sinqle  precision  routine)  and  USIMEC  which  inverts 
a  matrix  which  is  a  different  order  than  th<-  dimension  statement, 
are  also  used  throughout  KADILS.  Appropriate  system  calls  will  have 
to  be  applied. 

C.  Interference  Intensity  Calculation 

I.  By  substractinq  a  backqround  scattering  curve 

SI  (s)  •  (3)  -  B(s)  )/lf^ 

Background  line  is  calculated  from 

A.  Series  of  overlapping  exponential  curves 

1.  Starting  parameters  obtained  by  a  guess  at  the 
background  scattering  of  selected  values  of  s. 

2.  Starting  parameters  obtained  from  selected  values 
of  If^,  which  has  been  scaled  to  the  experimental 
total  intensity. 

B.  Summation  of  Gaussian  curves 

1.  Starting  p.irametors  are  determined  in  a  separate 
program  oy  a  trial  and  error  fit  to  the  If^. 

II.  By  sabtrac'ing  and  lividir.n  bv  a  b.ickaround  curve.  The  back¬ 
ground  is  c.alculated  is  in  i-'l  and  2  Only. 

si(8)  *  I^(s)  B(b)  ■*  1.0 

D.  A  few  Texas  Instrument  -  ASC  system  routines  are  used  for 
speeii  and  convenience.  These  system  routines  must  be  replaced. 
Enclosed  ifi  a  description  of  the  name  and  uses  of  the  ASC 
routines  used  to  facilitate  substitution. 
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»  1  ,  OijO 
p*&f  I 


sutPfiocp*^  lihhahv  seTEHnen 

03/ Jt /77 


SI»LL  ECo*M'*i'S  /  piTbI*  IKvESSIOk  /  C£TEfl>*JP*Ky 


I.  f"iMS 


►  'TRICES  st'SEC  et  CPLLP>>3I 

C‘lL 

si'-EC  ('.p.N.-.u.n 

PEtL*« 

*PCtP£*>T* 

C‘lL 

fct'L'P 

APCLfEMg 

C  *U 

CSI‘'EC  I  ) 

C*P 

*BCLiPE'‘T8 

c*ll 

cr'Txcd.e.^.x.o.n 

C  •  1  0 

iPcUPEMg 

►  'tBicf  s  sT*«Er  fit  »«. 

S  1 

c  *ll 

SI  'S  (i.p.K.x.r.T) 

Pf  «L  ‘P 

AfiCLT-E^TS 

c‘ll 

r  5  r  '  -  ( * ,  B.K  ,»•  ,r  ,  I ) 

•^E  'L  *0 

asclt'E*'  is 

c  'll 

C5I“Eh  (*,R,K,*>,r,I) 

c.e 

aPCLPE  P  78 

c*ll 

CC  SI*--  (i.P.N.-.r.D 

C  •  t  6 

'RCL*-!  78 

C»  lE 

w*LI7E''  Sj-i  Lt  >\£  ns  EC.*EIP*>Si 

hiTRIff'  «t»h^P  a,  c*LL►^S 

c  *LL 

LSI'-'C  .1  ,r.Tfo*5,  I  «.J«, 

IT  .  JT  ) 

hE  aL  •  p 

APG 

C'LL 

rLbl‘C(»,»,«,l  .r',!»6iS,!«.J«, 

IT.JT  ) 

P  f  A  L  •  P 

ARC 

C‘LL 

CLS!  C{«,v.-.i,r,:E6*S,I«.J», 

IT  ,  JT  ) 

C«-PLf  »  •» 

ARC 

C*LL 

crLS^r(i,v#»  f  r 

It,  JT  ) 

C "tPLE  »  •  1  t 

ARC 

‘■'tPlff*  «T»«»P  H»  E'«S 

C'LL 

L£I“EK(i,v,.,l ,r.lf 

IT.JT) 

WE  aL  ‘p 

•PC 

C'lL 

CLSI‘‘-(i,»,«,i  ,r.TE«'S.I«.J». 

IT.JT) 

w  E  * l  *6 

APb 

C'lL 

CLSI“«:(i,»,.,l  ,r,IEsiS,I*»J>, 

IT,  JT  ) 

C’pPlE  »  pp 

ARC 

C'LL 

CrLS-^(«(«,y,it,L,r,TES»S,I«.J«, 

IT,  JT  ) 

C'pPLf**it 

AEG 

?,  FLPPtSf 


T”?  S"iU  SivtBii  SETS  nr  «  I L  ^  f  L  5  ECL*Ti*kS.  iwtbT  « 
•-•TPn,  !•»>  P*Th.  It  ri‘  US*  (.SEC  fv*IL*Tl  * 

CETffi*-l»‘**'T, 

s.  »bce‘'Ekt  riE  it  iti»n^ 


I*>Pl  T 

i  ►iti.i,  "r  r'fTtlrlE*.*! 

f>  kiTuit  "r  r*^ST»sTS 

K  ‘.^‘’‘'E-  *r  b'.S 

„  t  a  •f  !*■  P  (“*T  PE  iEsP  ) 

I  PL't"  T«  fD*$»n.t  S’*fc|SE  *T  l.E*ST  <(K  Lf*'CTF, 

*ijEP(.  f 

t  «r(C  •■i’i.l* 
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»  <  .Ci,c 

I  aGI  t 


SLt>’acGH*>’  LiHR*wr  srifaneK 

OJ/l  I  /77 


St'^lL  t,  t-L  •  T  t  ('►  S  /  ‘*TPI*  iSvluSTf*'  /  CtUS^lK^sT  -  CCnTiNCEC 


H  7r  b?«T3 

C  Cf  7t»'“I*'*KT 

C*KFa*LWfn  bT“|iLT»K{«lj5  fCLATHNSl 

I  INPUT  “aTbI*  rf  C''EFFlClENTs  C  I  3  I  "sE  D  fI».JX) 

Y  YICT'K  »»«  *'»TcI»  NF  C»*S8T*NTS  Cl*'£T>3l8F-ECf  n  ,  JY) 

K  Kl.‘‘PfP  "F  H»«J  anO  C9is  'if  •«■  Tfl  be  L'SEC. 

L  rF  B(«"«  (I  SI-FC 1  9H  CruSlLSIBEB)  Bf'y  TP  i;SE, 


c 

C(  tEi.a'Ia.aa.T  (output  ) 

IFi.*S 

PL'C*  'f  ySaSaHiF  3TSBE*CE  fjo  Is 

LESC^H ) 

I  « 

K*-  iI-fsbIos  "r  •*• 

J  ‘ 

C*L  I^-fSBlos  Pr  ’»■ 

I' 

P*«  I“fSBloK  fY(  •»• 

J» 

f*l.  'I^fSBlo*.  f,-  'Y* 

t  I*-  JT- 

I'  *Nl. 

T  OA’l-I,  t*-vffcSi"s  •«  "tT{Hi,I».AST 

EvalUaT  i«s 

L'ESJKFC',  J'T  T«  I  ft- •  .  ►••F^Et,  'e*  “LST  STjLL  F£ 

FBf  S{  N  T  . 

FLIPPLC*-  I’  Is  •Nl  Y  A 
f.  '‘FPLEb 


U5l  5  I  ‘f  C  I 


*  f  ’  .?)  ... 

*(',*•) 

‘  (i,  '  1 

*(;>.;)  ... 

A (;,N ) 

•  •  • 

* (‘  .  '  ) 

*  (•  ,;) 

A  (S  ,S  ) 

•M ' .  n 

u  (  1  ,  o  ) 

►*  (?.;i 

PC?.*”) 

t  •  • 

'» ( s  .  1  ) 

P  (A  ,►  ) 

Call  si 

‘■I  c  fA.p.s.'-.r 

,  1 1 

’l'-fs 

•  • 

•  • 

.  • 

1  p  ( t . ' ) ) 

(»( 1 .?) 1 

1  - (  1  .s  )  1 

1  p  { ? . « )  1 

1  B  (  ^  (  2  )  1 

1  '•  {  2  ,  S  )  1 

A 

'*  '  . . .  1 

,  A  1  1  1  ...  1 

.  .  •  *  •  *  1  ...  1 

1  ■  (s  .  1  )  1 

t“(N.2) 1 

1  f  ( ‘  )  1 

•  • 

•  • 

•  • 

bFsTpIcTjnas 


T-(  »aTljc<S  ‘A-  ant  -f.*  aF(  CI!Tb«»iC  F»  FyECCtI'n  !»  T*-I 
5  pB-l.  T  INF  . 


>  t  .Uo'" 


SLEPReCfl*^'  LIH«*RY  NPTt«<HfK 

03/11/77 


SI*’lL  tCL*It''^S  /  TSvtBSI"^  /  Ct7EfiHlN*KT  •  Cf'NTiM.EC 


■*’,  "S".  •‘“'C  *0*  *Rf  Rt*L*P»  CC*‘PLE>«fl  OP 

C’^PLt»«lfc  IEPIkCisc  PS  the  'CAIL’  STME^'fM, 

.  *  i. ».  j  V  C 

I"t  -jtRlCtS  **•  aNT  "H*  “LST  he  P*C«tCI  Th*T  jj,  •4*  MST 
‘'ijST  PE  Ul'-Espi'^NFC  (k.N)  *sO  -e*  PLJ!  PE  T  T  **  E  •»  8  I  ^  ^  E  D  (S,p) 
•><  (m,k)  CEPE^Tt^G  •'<  ..HfThEB  TRE  p*TPIC£S  A^E  STORED  8» 
s-.S  *R  CUL-kS'.  f*A-PiE,  If  IT  IS  desired  TP  IKVERT 


4tl  I  3i 
I  ?  -El 


A*.c  ■*■  IS  !''rPP"FCTL»  nl^-E^SIO^Er  (3»3)»  T^fs  -a*  <«1lL 
C'nTaTn  ■'►•t  T  ‘•I  ►‘"ErSi 

3. •<.*.«.«.* 

asc  TRf.  51.  >^'ITt>.F  ’St^-FC*  •■UL  ATTE-PT  T"  I^'VEbT  Tre  ►aTBIi 


II  1 1 
FI 
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k' » I'.  I  1 


*•  !*  T  f  ti  I*  r » 

f'l/W/Tl 


•-irf*!*  ^'iiTlPLlC*Tlflk 


'  .  t  ►  ■  I  r  r  '  I  *■  T  ‘ 


C  ‘ll  »  1 1 .  J  .  A.M.C  ) 

c*ll  r*-^P''  (l.j.A.A.p.r ) 
C  ‘1.L  r-*  P*  ( I .  .  A  ,rl.c  ) 

C*LL  CC-“Pr(l.J.A,*.H.ci 


f.E  *L  •  U 
Bt  4  L  •  P 
ce^'PLf » •» 


jBCl^ekis 

iBCc^t  ^  T5 

*eCL''E  *.T! 


l.‘ Vf  MAC  I /E  r  ►•»TcI»-‘’Atwt»  “LLMPLICATICA  I 


call 

L“*'P'»  (1 

Hf  i  U  •  » 

ACG 

c  *LL 

r  L  »•  -  P  >  ( 1 

Bt  aL  • 

A  B  G 

c  ‘IL 

C  L  “  ^  p  *  ( 1 

C'yPi.E  1  .fl 
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V. 


GLOSSARY 


A(L)  -  Coordination  number  for  the  Lth  distance.  Corresponds 
to  in  Eq.  5. 

AA,  BB,  CC,  DD  -  parameter  for  calculation  of  x-ray  or  neutron 
atomic  form  factors.  Used  only  in  XRTI . 

AA,  BB,  CC  -  are  the  unit  cell  dimension  in  LDSIS.  CA,  CB,  CG 
are  the  direction  cosines. 

AK  (100,100)  -  is  the  coefficient  matrix  in  RArilLS.  After  solution 
of  the  problem  Ax  •  B  either  by  matrix  inversion  or  the 
e iqenvector-e igonva lue  method  AK  contains  the  inverse  matrix 
which  is  used  to  determine  standard  deviations  and  the 
correlation  matrix.  AK  also  serves  as  product  matrix  in 
IGN. 

AL(100,100)  -  scratch  array  use^l  to  store  coefficient  matrix  in 
INVERT. 

AM(100,100)  -  scratch  array  used  for  product  matrix  in  INVERT. 

A-MPMAX  -  IS  the  largest  value  of  the  amplitude  at  large  r. 

AP{3,1),  Br(3)  -  arrays  used  in  CURVE  to  obtain  initial  overlapping 
background  curve  parameters  from  the  Tf^. 

APACK  ( 1 0 , 000)  ,  BRACK  ( 1  00)  -  paclied  vectors  containing  AK  and  BK 
in  INVERT. 

AT(100,100)  -  is  the  transpose  of  the  eigenvector  matrix  CK. 

AZ<M,NC),  BZ(M,NC),  CZ(M,NC),  DZ(M,NC)  are  the  exponential 

shape  parameters  of  the  mth  se<Tment  of  the  NCth  background 
line . 

B(L)  -  the  B  values  of  the  Lth  distance.  Corre'^ponds  to 

in  to.  5. 

BAR(K,NL)  -  is  a  normalized  scattering  factor  for  atom  NL  at  the 
Kth  value  of  S.  It  can  be  expressed  as 

f,.j  (s)  ^vr”  f^ 
uc 
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BK(IOO)  -  IS  the  Weighted  residual  vector  in  RADILS.  After  matrix 
inversion  or  the  eigenvector-eigenvalue  method  the  parameter 
changes  are  storru  in  BK. 

Bt.(100)  -  scratch  array  used  to  store  residual  vector  BK  in  INVERT. 

RLK(-IOOO)  -  is  a  scratch  array  used  in  ICN  and  ultimately  contains 
the  linear  parameter  changes. 

BROADN(600)  -  scratch  .irray  used  to  calculate  Iroadened  total 
intensity  of  a  crystalline  topology  in  BROAD. 

CDIS (  3500 , 10)  -  IS  the  coordination  number  for  each  type  of  distance 
in  bins  of  0.01  A  out  to  DISM.X. 

CH  -  the  camera  heicjht  for  the  current  data  set.  L'sed  only  in  ED  I L . 

CK(100,100)  -  Coefficient  matrix  AK  is  stored  in  CK  in  IGN.  After 
call  to  EIGRS  the  eigenvector  matrix  ,X  is  stored  in  CK . 

CO,  FA,  .SK5  -  parameters  used  in  the  aaussian  summat  ion  background  line 

2 

Bis)  •  SK5  (  5  A  expl-H  s  )  +  FA.s  ♦  CO) 

'  n  ^  n 

m*  4 

CSIS(N)  -  the  combineti  reduced  intensity  corres[>onding  to  si(s)  in 
EJ.  5  at  slN).  For  multiple  data  ranges  the  individual 
intensities  are  combined  in  SCBSIS. 

DAMP  -  Damping  factor  applied  to  DRIS  in  SUBSIS,  and  used  in  the 
refinement.  DAMF  corresponds  to  i  in  FQ. 

2 

DA>tA  -  the  value  Of  a  in  FQ.  7  can  be  set  equal  to  -DAMyX.r  for  the  final 

radial  d  i  .s  t  r  i  bu  t  ion  plot.  This  variation  of  i  with  respect  to  r  was 
introduced  in  order  tc'  diminish  the  effect  of  random  errors  at 
larger  s  on  the  features  of  the  RDF  at  large  r,  which  arise 
primarily  from  the  low  angle  data. 

DELIL)  -  the  Lth  observer!  value 

Ir.  (r^)  for  1<I.  N0BR(1) 

r(  (r  )  -  .  )  for  NOBS  ( 1 )  ♦  1<L<N0BS  ( 1 )  ♦NOB.S  ( 2  ) 

L  *  o' 

._Sl  (s)  NOBS  (1)  ♦NOBS  (2)  +1<I.<N0BS  ( 1 )  ♦NOBS  ( 2  )  iNOBS  (  3 1 

/■.Si(s)  NOBS(l)^NOBS(2)^NOBS(3)^l  L - 

etc . 
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DET  -  the  value  of  the  determinant  obtained  from  the  ASC  system 
matrix  inversion  routine. 


Ofliw  -  IS  the  total  standard  deviation  after  the  current  refinement 

/NKHO - - 

OKUW  -  I  ^  DEL(  I  ) '^/ (NRHO-NPARR  +  NCOS 

DEX(20)  -  array  used  in  Cl’RVE  to  obtain  initial  value  of  c  in  expia  ♦ 
from  a  3-point  fit  to  initial  background  points. 

DIJd.ICT)  -  derivative  matrix  DIJ(200,50).  It  represents  the 

derivative-  of  the  ICT  parameter  evaluod  at  the  Lt h  observed 
%’a  lue . 


OIND(200,S0)  -  IS  the  derivative  matrix  ot  independent  parameters 
when  constraints  are  used  on  the  structural  paramr-ters. 


DISMN 


DISMX 


is  the  minimum  distance  at  which  LDSIS  calculates  the 
longer  distance  il.D)  contribution  in  I.bnI.s.  .<0  refers  to 
short  ordered  distances. 


8i'(s)  -  31(S) 


(si 


SD 


•NDIST 

0 


sij^D<s' 


DISMX 


) 


DISMN 


t  h<’ 


IS  the  maximum  v.ilue  of  r  where  the  long  distance  contributors 
are  calculated  in  LDSIS. 


DK(IOO)  -  Weighted  residual  vector  BK  is  stored  in  DK  in  ICN . 

DOBS(K)  -  Increment  of  N  for  the  Kth  observation  ranocs.  Read 
in  only  for  K  «  1  4  2. 

DS  -  the  increment  in  As 


DSIS(N)  “  CSIS(S)  -  TSIS(N).  This  corresponds  to  si'(s)  in  EO.  5. 
DlMdOO)  -  scratch  array  used  m  USIMEC,  ASC  system  matrix  inversion. 
DUMl  -  not  used. 

DX ( 1 2 ) ,  DY ( 1 2 ) ,  DZ(12)  -cell  translation  for  x,  y  and  z  respectively 
on  symmetry  transformat  ion  card  in  LDSIS. 
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ED(27) 


EPS 


-  the  array  used  to  print  out  the  final  normalized  values 
of  the  data  in  SUBTR. 

1.0  -  1.5  X  r/D  +  0.5  (r/D)^  r  ^  n. 

EPS  is  a  spherical  particle  size  function  which  corrects  the 
relative  frequency  of  occurrence  of  an  interatomic  vector  in 
a  structure  of  unbound  dimension  to  that  which  pertains 
within  a  sphere  of  diameter  D. 


ESIS(600)  -  scratch  matrix  used  to  smooth  DSIS  in  SUBSIS. 

EV(IOO)  -  IS  initially  the  vector  of  eiqenvalues  in  IGN.  Subsequent 
to  call  to  EIGRS  the  inverse  oifjenvaluer  are  stored  in  EV. 

EV.MIS  -  minimum  eiqenvalue  accepted  in  filterinq  process  in  IGN. 

EVSAV(IOO)  -  the  eiqenvalue.s  in  IGN  are  stored  in  EVSAV. 

FA  -  Parameter  used  in  rampinq  function  RAMPIS).  If  IDX(14)<0 
RA.yp(s)  -  1.0  -  EXr(fA.5)  otherwise  RAMP{S)  -  FA.s.  If 
IPX  (20)  <0  RAMr(s)  -  -100  1  This  function  could  be 

experimented  with,  to  obtain  the  best  shape  for  backqround 
f  1 1 1 1  nq . 


FACTOR  -  not  used. 

FBAR(N,NL)  -  IS  the  scatteriiiq  factor  coefficient  for  the  NLth 

atom  pair  at  the  Nth  value  of  s.  It  can  be  expressed  as 


f  (s)  •  f  (s)  •  cos  -  n.'s) 

A  B  BA 


the  f  and  h  '•epresent  the  maqnitude  and  phase  of  the  coherent 
atomic  scattorinq  factor.  The  cos  term  is  only  used  in  ED. 


IS  the  sharpeninq  function  I  (s) .  It  is  calculated  for  XP 

s 

and  ND  in  XRTI  from  the  scatter inq  factors  and  it  is  read 

from  file  for  ED.  I  (s)  »  1  where  uc  represent  the 

.  .  f  ...  ®  uc 

unit  of  composition. 


FC  (N) 


FF(I) 


-  the  value  of  the  Ith  point  of  the  initial  backqround 
for  the  current  data  ranqe.Useii  in  CURVK  to  generate 
initial  background  parameters. 

FJ  -  DS/(2n^) 

FOBS(K)  -  the  first  r-value  for  the  Kth  observation  range. 

FOBS(K)  is  different  from  zero  only  for  K  »  1  (inner 
region  of  the  RDF)  and  for  K  ■  2  (outer  reuion  of  the  RDF). 


FT(N,NC)  -  Continuous  compiosite  background  for  the  NCth  data  set 

of  the  Nth  r-value  made  up  of  weighted  overlapping  exjxinential 
background  segments. 


FT(NC)  -I]w,^(s)  exp  (a. 

^  ^Tn  ^ 

m 


km 


GK(100,100)  -  scratch  array  used  in  IGN  for  the  eigenvectors  and 
dimensioned  in  RADIL.S. 


lA(lOO)  -  scratch  array  used  in  USIMKC  the  ASC  system  matrix  inversion. 

ICT  -  is  an  index  of  the  parameters  in  the  derivative  matrix  DU. 

IDX(I)  -  the  flag  array  used  to  control  the  various  options. 

IN(I,NC)  -  the  integer  corresponding  to  SNN(I,NC). 

T  PnwvR 

iroWF.R  -  the  power  of  N  in  the  RDF  plot.  ^ .  j.  '  • ^ 

IS  plotted.  The  default  value  is  1. 

ITYPF.  -  Flag  for  treatment  of  the  background.  If  IDX(18)  >  0  the 

empirical  background  is  subtracted  to  obtain  SIS.  If  IDX(18)<0 
SIS  is  obtained  by  suntract inu  and  dividina  bv  the  background. 

KI(N)  -  the  variation  flag  for  the  Nth  parameter.  Refer  to  refinement 
section  for  clarification. 


KRAVG  -  Running  average  parameter.  If  greater  than  0,  a  (2*  KRAVG  ♦  1) 
point  running  average  on  the  t^otal  intensity  is  performed  so 
that  the  correlation  length  is  KRAVG»DS.  If  IDX(24)<0  DSIS  is 
smooth  in  three  stages.  The  data  from  0.05  -  5.0  are  not 
altered.  On  the  first  smoothing  cycle  ESIS  is  smooth  with 
"IDX(23)"  passes  of  a  (2*KRAVG  1)  fxiint  square  wave  and  the 
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data  in  DSIS  from  5.0  to  N’SIS2  are  replaced.  On  the  second 
smoothinq  cycle  the  already  smooth  ESIS  is  aqain  smooth  as 
before  and  the  data  in  DSIS  from  15.05  to  NSIS2  are  replaced. 

Ml  -  set  to  100  -  not  used  in  proqram. 

HAXPAR  -  maximum  nu.mber  of  linear  parameters  (eiqenvalues)  allowed  in  ION. 

MTITL  -  Title  card  for  job. 

Nl  -  the  number  of  variables  +  backqround  I.aqranqian  multipliers. 

It  IS  the  order  of  the  matrix  AK,  CK,  AT,  X. 

NADP  -  Klaq  for  reiut roduc inq  the  short  distance  peaks  (without 

termin.ition  errors  or  errors  due  to  non-constant  coefficients) 
into  the  plot  of  the  RDF.  S'ADD  -  0  distances  are  not  re¬ 
introduced.  NADD  ■  1  distances  are. 

NA.ME2,  .NA-MF-l  -  is  the  date  of  the  job  which  created  the  ED  scattering 
factor  file. 

SAMEl  -  also  represents  the  JOP  ID  of  the  job  which  created  the  ED 
scattering  factor  file. 

SA.yF2  -  dash(-) 

NA,MKl,  NA.HE3  -  alphanumeric  atom  name  for  ED  scattering  factor 
coef  f icients. 

NATAS  -  number  of  atoms  in  the  as^mimetric  unit  in  LDSIS. 

NATTY  -  number  of  distinct  atom  types  in  sample  in  LDSIS. 

NC  -  the  index  which  normally  corresponds  to  the  data 
range  currently  being  processed. 

NCAS  -  the  number  of  ranges  of  data  to  be  processed.  The  maximum  number  Of 
ranges  is  -J .  The  present  version  of  the  program  has  been  tested  for 
NCAS  »  1  only. 

SCIN  -  flag  for  the  type  of  background  input.  If  IDK  (17)<0,  NCIN  »  1 
and  AZ’s,  BZ's,  CZ's  and  DZ's  for  the  exponentials  are  read  in. 

If  IDX(17)  ^  O.NCIN  «  0  and  AZ’s,  BZ’s,  CZ's  and  DZ's  are 
calculated  from  points  on  the  initial  backqround  which  are 
read  in. 
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NCOEKF 


-  the  number  of  distinct  coefficient  types  to  be  read  in. 
Used  in  EDIL. 

NCO.V  -  number  of  constraints  imfwsed  on  the  structural  parameters 
in  CON. 

NCV(NC)  -  The  number  of  backqround  curves  making  up  the  background 
for  the  NCth  d.it.i  r.in<?e.  JK  is  often  used  to  be  the  number 
of  curves  for  the  current  data  r-inue. 

NDA  -  the  number  of  distinct  atom  types  in  the  unit  of  composition. 

NDIST  -  the  number  of  short  distances  for  which  parameters  will 
be  read  in. 

NKMT  -  the  format  to  be  used  for  the  input  of  experimental  x-ray 
or  neutron  intensities. 

NG  -  number  of  gaussians  used  in  Gaussian  summation  background  line 

ng 

B(s)  -  Sk5  (  ic 
n*  1 

NIV  -  is  the  number  of  independent  variables  when  the  structural 
parameters  are  constrained. 

SNTYPECIOO)  -  atom  type  identification  for  atoms  in  as^irenetric 
unit  in  I.DS  I S  . 

NOn.S(K)  -  the  number  of  observations  in  the  Kth  observation  range. 

K  “  1  (inner  reuion  of  '  <  t )  * 

K  •  2  (outer  region  of  G’(r)) 

K  -  3,'l,5  (overlap  of  various  data  sets.) 
only  NOBS ( 1 )  and  NOBS (2)  are  calculated. 

Nr(L)  -  the  pair  type  of  the  Lth  distance. 

NPAR  -  number  of  flags  on  KI  card. 

NPARM  -  location  of  flag  on  KI  card  indicating  background  treatment 
-  the  number  of  variables  being  refined. 
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NPARR 


NRHO  -  Total  number  of  observations  -NOBS(N).  The  derivative 
matrix  is  NRHO  X  NPARR. 

NSD  -  number  of  ranges  over  which  a  standard  deviation  is 

calculated  NSD  “  2  (NCAS  -1).  2  is  for  the  inner  and 

outer  region  of  the  RDF  and  (NCAS-1)  is  the  number  of  overlap 
overlap  regions  in  the  si(s)  curve  for  NCAS  ranges  of  data. 

NSISl  -  the  index  of  the  largest  experimental  s-value. 

NSIS2  -  Index  of  the  largest  s-value  used  in  the  least-square 
ref inement . 


NSIS2T  -  number  of  ED  scattering  factor  coefficient  points  on 
input  file  tc  EDIL. 

NSUB  -  Flag  for  displaying  the  intensities  and  RDF  before  the 

current  refinement  cycle  begins.  NSUB  “  0  for  no  display 
and  NSUB  •  1  for  a  display. 

NSYM  -  nu.mber  of  s>Timetry  transformations  in  LDSIS. 

NTITL(I)  -  Title  array  describing  the  current  data  set  used  in  EDIL. 


NXl  -  10  -  DISMX/AA 
NX2  -  10  DISMX/AA 


define  the  number  of  unit  cells  in  the  x 


direction  in  which  to  search  for  contributing  atoms  in  LDSIS. 


NY  I  -  10  -  DISMX/BB 

NY2  ■  10  ♦  DISMX  BB  ”  define  the  number  of  unit  cells  in  the  y  direction 
in  which  to  search  for  contributing  atoms  in  LDSIS. 


NZl  -  10  -  DISMX/CC 
NZ2  -  10  ♦  DISMX/CC 


define  number  of  unit  cells  in  the  z  direction 


to  search  for  contributing  atoms  in  LDSIS. 

PCTSG  -  A  constant  error  expressed  as  a  percent  of  the  sharpening 

function.  Used  in  EDIL  to  generate  an  SG  array.  PCTSG  *  0.005 
corresponds  to  0.5»  of  the  sharpening  function. 

PFR(IOO),  PEBSAV(IOO)  -  are  the  fractional  linear  residual  reduction. 
PLOTA  -  ASC  system  scratch  array  used  for  plotting. 
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PSZ 


-  is  the  radius  of  a  spherical  ordered  region  with  a  particular 
bonding  topxslogy. 

PTDS  -  the  plate  to  sector  separation  for  the  current  data  set. 

Used  only  in  EDIL. 

Q(100,100)  is  a  matrix  (NPARR  X  SPARER)  which  is  composed  of  the 

coefficients  of  the  dependent  parameters  in  the  linear  constraint 
relations.  The  NIV  independent  parameters  occupy  the  first  NIV 
rows  of  Q  with  ones  in  the  appropriate  columns. 

OISV ( 100, 100)  -  inverse  of  Q  equal  to  J  matrix  in  constraint  writeup. 

QT(100,100)  -  transfKJse  of  0  Described  as  in  constraint  writeup 

and  relates  the  dependent  and  independent  derivatives. 

RA.Mr(S)  -  A  ramping  function  used  to  make  the  observed  background 
more  amenable  to  fitting  with  a  sum  of  exponent la 1 s . 

R.\Mr(S)  -  1.0  -  exp  (FA*S)  if  IDX(14)  <  0 

or 

-  FA*s  if  IDX(14)  >  0 

RDD,  WDD,  RDDZ,  VTDDZ,  A-MPMAX  -  An  expKDnential  fit  of  the  function 

A.Mr  ■  AMPMAX  (1.0  -  EXP  (Ax  ♦  Bx  x  r)  to  these  parameters  is 

used  in  LDSIS  to  determine  the  thermal  parameters  for  the 

2  2 

longer  distances  in  the  intensity  curve  EXP  (-s  AMP  /2.0). 

RDD  and  WDD  are  the  distance  and  amplitude  for  the  first  bonded 
distance. 

RDD2  and  WDD2  -  are  the  distance  and  amplitude  for  any  other  distance 
for  which  a  good  approximation  for  the  amplitude  is  known. 

RMAX  -  the  maximum  value  of  r  (in  A)  to  be  plotted  in  the  RDF. 

Used  only  in  SUBTR. 

RSMIS  -  minimum  residual  reduction  allowed  in  filtering  process  in 
IGN  for  a  particular  eiqenshift. 

RZ(L)  -  the  r^  of  tt.e  Lth  distance,  corresponding  to  r^^  in  Eq.  5 

s(N)  -  the  Nth  s-value  s(N)  «  8(1)  *  N*DS  where  s(l)  is  the  initial 
value  of  s  and  DS  the  s  increment. 
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SCTI  -  Scale  factor  tor  total  intensity  and  standard  deviation. 
Can  be  used  to  scale  date  prior  to  analysis. 


SD(IDS)  -  13  the  value  of  the  standard  deviation  after  the 
current  refinement  cycle  for  the  IDSth  range  of  DKL. 


SD(IDS) 


■r 


NOBS (IDS) 


E 

L*1 


DEL (L)  /NOBS (IDS) 


SDOLD(IDS)  -  IS  the  standard  deviation  for  the  IDSth  region  of 

the  least  squares  minimization  before  the  current  refinement 
eye 1 e . 


SDOLD(IDS) 


V 


NOBS (IDS) 

E 

L“  1 


DEL (L)  /NOBS (IDS) 


SE(NC)  -  the  scale  factor  (SK/SK3(NC))  necessary  to  bring  the 
the'  retical  atomic  intensity  to  the  level  of  the  NCth 

data  set  fur  the  subtraction  in  EDI L .  If  SE(NC)  ■  0.0, 

no  atomic  theory  is  subtracted. 

SERST(NC)  -  the  minimum  s-value  for  the  NCth  data  ranue. 

STT  -  the  fraction  of  the  shift  vector  calculated  in  the  least- 

squares  refinement  which  will  actually  be  applied.  SET  •  1.0 

jiv«*s  the  full  shift  SET  •  0.5  gives  half  of  the  shift,  etc. 

.sC.(N)  -  The  statistical  un./er  t  a  i  n  t  y  of  TI(N,NC)  at  s(N).  Eor  XD 
and  .ND  it  i:.  th*'  count  in-i  statistics.  Eor  ED  it  equals 
a  fM’rcent  uncc'rtainty  in  TI(N,NC). 

SH(IOO)  -  contains  accumulated  sum  of  parameter  shifts. 

SHIET  -  maximum  linear  parameter  shift  (oigenshift)  allowed  in  filterii 
1  rijc«*ss  in  ir.N. 

SIS(S',NC)  -  the  reduced  intensity  for  the  NCth  data  set  at  s(N)  . 

.‘‘•K  -  An  overall  scale  factor  used  to  bring  the  theory  to  the  same 
scale  as  the  experiment.  SK  corresponds  to  SK  in  EO.  7. 


SK3(NC)  -  the  scale  factor  for  the  NCth  -i.itn  set  used  to  brlnq 
the  various  data  sets  to  the  same  level.  Used  only  for 
multiple  sets  of  overlapping}  data  or  to  scale  experimental 
data  to  yield  an  SK  a  1 . 

SK5(NC)  -  the  overall  scale  factor  of  the  background  for  the  NCth 
data  set. 

SKK  -  A  scale  factor  lor  the  total  intensity  plot  in  SUBTR.  If 
SKR  '  0.0,  TI(N,NC)  and  SKR'TUN.NC)  are  both  plotted. 

SLST(NC)  -  the  maximum  s-value  for  the  NCth  data  set. 

SNN‘(I,NC)  -  the  midpoint  for  the  Ith  exponential  segment  of  the 
NCth  dat.i  set.  The  first  segment  of  the  NCth  data  set 
ranges  from  SFRST(NC)  to  .SNN(2,NC),  the  last  segment  of  the 
NCth  data  s.-t  ranges  from  .SNN(NCV(NC)  -1,  NC)  to  SLST  (NC)  , 
and  the  .Mth  segment  {yy  .^jcv  (NC)  ) from  SNN(y-l,NC)  to 
SNN (Mtl ,NC) . 

SOL(N)  -  are  the  dependent  parameter  changes  obtained  from  the 
independent  parameter  changes  when  constraints  on  the 
structural  parameters  are  employed. 

SS ( I )  -  the  s-values  for  the  Ith  point  of  the  initial  background 

for  the  current  >!ata  range  used  in  Cl.' PVK  to  generate  initial 
background  parameters. 

Symmetry  transformation  matrix  and  translations  vector  for  the  Mth 
symmetry  operation  on  the  LTth  atom  in  the  asymmetric  unit 
yielding  transformed  coordinates  XAA,  YAA,  7AA. 


FXl (M) 

FX2 (M) 

FX3 (M) 

XX{I.T) 

DX  (M) 

XAA 

FYl (M) 

FY2 (M) 

FY3 (M) 

YY (LT) 

•f 

DY  (M) 

« 

YAA 

FZl (M) 

F7,2  (M) 

FZ3 (M) 

ZZ (LT) 

DZ  (M) 

ZAA 
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TI(N,NC)  -  the  intensity  for  the  Mi'th  Jata  ranqe  at  the  Nth  s-value. 

TIA(N)  -  the  total  theoretical  atomic  intensity  at  the  Nth  s-value. 

TSDOLD  -  IS  the  total  standard  deviation  over  the  NHRO  points  of 
DEL  before  the  current  refinement  cycle. 

VNhkr-  ~ 

V  DEL(L) ^/(NHRO-NPARR) 

L-1 

TSIS(N)  -  the  theoretical  reduced  intensity  at  »(N). 

TSISOIS)  -  the  theoretical  reduced  intensity  at  «(N)  corresponding 

to  the  long  distances  generated  in  l.DS  I S .  This  array  is  zeroed 
if  no  lung  distance  contributions  are  calculated. 

U(3),  V(3),  W(3)  -  are  arrays  used  in  Cl'RVF  to  indicate  the  s  value, 
observed  background  and  calculated  backqround  obtained  from  a 
3-point  fit  to  initial  background  points. 

V(IOO)  -  scratch  array  used  in  ION  to  calculate  the  product  of  the 
transposed  eigenvector  matrix  and  residual  vector. 

WAVENO  -  (2k)  or  (4n/X)  for  the  current  lata  set.  Used  only  in  EDIL. 

WOBS(K)  -  the  weight  of  the  Kth  observation  range. 

WT  -  the  weighting  function,  which  is  used  ti-  fcrr  the  .'(Mnbined  intensity 

(CSIS'  from  the  individual  intensities  (SIS's)  when  multiple  ranges  are 
used.  WT  IS  0.5  in  an  overlap  and  1.0  otherwise.  This  combination 
is  performed  In  SUBSIS.  Used  only  wnen  NCA.s  e  1. 

WT(300)  -  is  the  weighting  applied  to  Insure  proper  coordination 
number  calculation  in  LOSIS. 
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WTNG(IOO)  -  IS  a  weight  vector  generally  equal  to  the  reciprocal 
of  the  square  root  of  the  diagonal  coefficient  matrix 
elements.  It  is  used  to  normalize  the  coefficient  matrix 
diagonal  elements  to  unity. 

WT1,WT2  -  The  weighting  functions  used  to  combine  the  exponential 
segments  lo  obtain  the  background  for  the  NCth  data  set. 

WTl  »  WT2  *0.5  foi  overlapping  segments  and  WTl  *  1.0  in 
the  non-overlap  regions.  WTl  and  WT2  correspond  to 
in  EQN.  8.  This  combination  is  performed  in  SUBBG  and 
SUBS  IS. 

.X(N1,N1)  -  eigenvector  matrix  is  obtained  from  EIGRS  in  this  array, 

XCHG{I)  -  is  the  change  calculated  from  the  current  refinement 
cycle  for  the  Ith  parameter. 

XLAB  -  the  amplitude  of  vibration  for  the  current  short  distance. 
The  B-value  is  calculated  by  B(L)  »  -XLAB  x  XLAB/2.0. 

XNEW(n  -  18  the  new  value  of  the  parameter  1  after  the  current  cycle. 

XOLD(I)  -  is  the  value  of  parameter  I  before  the  current  refinement 
cycle . 

XSD(I)  -  is  the  standard  derivation  for  the  Ith  parameter 

XSD ( I )  -  ^AK(I,I) ‘DEUW 
where  AK  is  the  inverse  matrix. 

XX(I,J)  and  YY(I,J)  are  the  scratch  arrays  used  in  the  call  to  the 

plotting  routine  CCPLT  to  represent  the  abscissa  and  ordinate 
of  the  Ith  point  of  the  Jtn  curve. 

XX(300)  YY(300)  ZZ(300)-  are  the  fractional  coordinates  of  the  atoms 
in  the  asymmetric  unit  in  LDSIS. 

YSD(N)  -  are  the  dependent  parameter  standard  derivations  obtained 

from  the  independent  parameter  deviations  by  the  transformation 
matrix  QINV. 

-  is  a  rejection  diagonal  matrix.  It  has  a  value  of  1  for  an 
eigenvalue  which  will  be  accepted  and  zero  otherwise. 
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IV.  PLOTTI.NG  PACKAGE  -  PLOTT 
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A.  Purpose 

The  purpose  of  this  subroutine  is  to  obtain  directly  on  a  line 
printer  a  crude  plot.  The  subroutine  will  plot  any  qiven  number  of 
arrays  on  a  variable  size  surface  limited  only  by  the  printer  capacity. 

The  subroutine  does  not  disturb  the  .arrays  being  plotted.  Alternate 
entry  points  allow  for  plotting  every  2nd,  3rd  or  5th  point,  or  a 
calcomp  plot. 

B.  Usage  of  the  Program 

1.  Operational  Procedure:  Usable  as  a  subroutine. 

2.  Input  data:  The  calling  8c<}uence  is: 

SUBROUTINE  PLOTT  ( ,XX  ,  YY  ,  NDATA  ,  NDMAX  ,  I  SYMBOL ,  NF , 

XLI NE  ,  M.\  ,  YLINE  ,MY  ,  SLS  ,  NCL ,  M.M ,  LI, ,  AREA ,  YSCALE)  . 

where , 

XX  (J, I)  and  YY(J,I)  arc  the  coordinate  vectors 
of  the  j-th  point  of  the  I-th  array  to  be 
plotted.  (e.q.  1=1,  to  plot  one  array). 

NDATA(J)  IS  the  number  of  points  on  the  J-th  array. 
NDMAX  IS  the  maximum  number  of  points  in  the 
largest  array  to  bo  plotted. 

In  the  dimension  statement: 

X.X(J,I)  and  YY(J,I):  J  must  be  the  value  NDMAX. 
(c.g.  NDMAX  =  100,  sec  NDATA(J)) 

ISYM{J)  is  the  s>T7ibol  used  in  the  J-th  array.  The 
ISY.MBL  array  may  be  defined  by  a  Hollerith 
statement  (e.q.  ISYM  (1)  » 

NF  is  the  number  of  arrays  to  be  plotted.  In  the 
dimension  statement: 

XX (J, I),  YY(J,I),  NDATA(I),  and  ISYMBL(T): 

I  should  be  equal  to  or  greater  than  the  value 
NF ,  (e.q.  NF  =  1,  to  plot  one  array) 

XLINF.(I)  is  the  value  for  the  1-th  x-reference  line, 
(e.q.  leave  undefined). 


MX  is  the  number  of  x-reference  lines.  In  the 
dimension  statement: 

XLINE(I);  I  should  be  equal  to  or  greater 
than  the  value  MX.  (e.g.  MX  =  0,  if  no 
x-reference  lines  are  desired) . 

YLINE(l)  is  the  value  for  the  I-th  y-roference 
line.  (e.g.  leave  undefined). 

MY  IS  the  number  of  y-reference  lines.  In  the 
dimensions  statement: 

YLINEd)  :  I  should  be  equal  to  or  greater 
than  the  value  MY.  (e.g.  MY  =  0,  if  no  y- 
reference  lines  are  desired). 

NLS  is  the  number  of  lines  used  in  the  plot.  In 
the  dimension  statement: 

AREA(30,J)  and  ySCALF.(J):  J  should  be  equal 
to  or  greater  than  the  value  NLS.  (e.g. 

NLS  ■  56,  for  one  page).  If  SLS<0,  the  tic 
marks  of  the  calcomp  plot  will  not  be  labeled. 
The  x-dimension  of  a  calcomp  plot  is  NLS, '10 
inches . 

M.M  control  integer:  MM  ■  1  plots  across  the  page: 

MM  ■  2  plots  down  the  p.tge. 

LL  control  integer:  LL  »  1  linear,  LL  «  2  semi-loo, 
LL  »  3  square  root . 

ARKA(30,J)  IB  the  vector  which  spans  the  plot 

surface.  This  vector  is  in  the  call  list  to 
satisfy  the  variable  dimensioning  requirements. 
YSCAL(J)  IS  the  scale  scalar  running  down  the  page. 
This  scalar  is  in  the  call  list  to  satisfy 
the  variable  dimensioning  requirements. 

3.  Output  Format;  Not  necessary  for  user  to  provide  output  format. 

*.  Caution  to  Useri  The  number  of  lines  (NLS)  drastically  alters 

the  storage  needed  m  plotting  because  of  the  array 
AR>:a ( 30 ,NI.S)  .  NLS  >  56  usually  is  near  to  exhausting 
a  page  of  output.  In  cases  where  the  user  wishes 
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20  or  more  pages  in  one  plot  the  dimension  of 
AREA  will  be  significant. 

The  user  must  alter  the  restrictive  statement 
in  PLOTT  if  NLS  is  needed  to  be  greater  than  512, 
if  MX  is  needed  to  be  greater  than  32,  or  if  NF 
is  needed  to  be  greater  than  16.  The  restrictions 
are  only  intended  to  prevent  unknowledgeable 
programmer  errors. 

5.  Error  Provision:  Error  stops  are  printed  in  output. 

Most  previously  used  problems  have  been  that 
the  values  for  variables  in  the  CALL  list  are  not 
what  the  programmer  intends.  Print  out  all  variables 
in  the  CALL  list  as  a  first  step  in  seeking  errors 
and  check  list  against  SL'BROCTINE  list. 

6.  Alternate  entry  points:  All  entry  points  have  the  same  parameter 

list. 

PLOTT2  -  will  plot  every  2nd  point  of  each  vector. 

PLOTT3  -  will  plot  every  3rd  point  of  each  vector. 

rLOTT5  -  will  plot  every  5th  of  each  vector. 

CCPLT  -  will  also  yield  a  calcomp  plot  of  dimensions 

(NLS/10)  'by  10"  and  by  making  NLS  and  NCL 
negative  the  labels  of  the  axis  tie  marks 
are  omitted. 

7.  Subprograms  Required:  Not  applicable. 

8.  Accuracy ;  The  plotting  accuracy  ir  limited  by  the  discreteness 
of  the  print  spacing  or  the  precision  of  the  calcomp  plotter. 

C.  Mathematical  Method 

The  maximum  and  minimum  values  of  each  coordinate  are  found.  The 
data  arc  than  transformed  by  an  appropriate  transformation  such  that 
the  values  range  from  0  to  "NCL"  for  the  horizontal  coordinate  and  from 
0  to  "SI.S"  for  the  vertical  coordinate.  These  are  rounded  and  truncated 
into  integer  values.  The  integer  values  are  then  used  to  locate  the 
Hollerith  symbol  in  memory.  If  this  location  is  already  occupied  by  a 
symbol,  an  equal  sign  is  substituted  for  the  s^mibol.  The  entire 
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Hollerith  array  than  represents  the  plots  and  is  consequently  printed. 

For  calcomp  plotting,  it  is  assumed  that  the  arrays  to  be  plotted  will 
be  ordered  into  either  increasing  or  decreasing  values  of  the  abscissa 

D.  Kxa.-nple:  For  a  15A~^  plot  at  1a"VinCH,  NLS  -  150  XLINE(l)  -  0, 
XLINE(2)  -  15.0 

The  program  will  put  out  both  a  line  printer  and  calcomp  plot. 
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VII. 


LKASl-SQUARK.S  MKTHOD 


SYMIKJLS  : 

V  -  function  to  be  minimized. 

K  -  a  function  of  the  x  dependent  variables  describini^  the 
observat lonal  parameters . 

-  observed  value  of  K. 

K  -  calculated  value  of  K. 

-  the  residual. 

D  -  the  derivative  matrix  •  3K./Dx^  for  the  ith  observation 
and  the  jth  ;'aramet*'r. 

Vi ^  -  wei.ght  for  the  ith  observation. 

A  -  tht‘  coefficient  matrix  D  P. 

■\\  -  the  shitt  vector  . 

ux’  -  normalized  par, meter  chan.je. 

Ax"  -  normalized  eigenshifts  (linear  combination  parameter  changes). 
B  -  residual  vector  wei<)hted  by  derivatives  D^F. 

N  -  norma  1 1 zat ion  weight  matrix 
T  -  eiuen vector  matrix, 
t  -  transpose  operation. 

•  -  eigenvalue. 

X  -  dependent  parameters. 

V  -  independent  parameters. 

z  -  filtering  diagonal  matrix  of  ones  and  zeros. 

-  total  standard  deviation, 
s  -  standard  deviation  for  parameters. 

S  -  correlation  matrix. 


E 


reduced  independent  derivative  matrix 


E 


Ik 


m  -  number  of  constraint 
M  -  number  of  observations. 

n  -  number  of  unconstrained  dependent  variables, 
k  -  number  of  independent  variables. 
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Matrix  Equations 


V 


M 

N' 

1-1 


“V< 


We  wish  to  ir.inimiz«^  this  equation 


3V 
3>! . 


31-'^ 

0  =  E  K  (F°  -  y^) 

“  1  1  i'  DX^ 


M 

V 

1-1 


( 3-1 , 2 . . . n) 


This  qives  the  set  of  n  normal  equations  in  n  unknowns. 


When  K*"  is  not  linear  it  must  be  approximated  by  a  Taylor  series 


in  where  hiqher  order  terms  are  neqlected 
1  '^2 


:)F(aj  .a^.  . 

K(x, ,x,. . .x^)  -  F(a, ,a,. . .a  )  ♦  - 


x 


1 


31  la^/a2f  *  • 


Jx 


ix. 


When  the  ^  '  s  good  approx  imaf  ions ,  application  of  the  least 


j 

squares  will  yield  values  for  the  parameter  changes  '.^y  where 


X  -  X  1  j  ♦  X 

inew  3 old  3 


In  matrix  notation  the  normal  equations  can  be  written  as 


Aax  -  B 


where 


A  ■  D*'D  B  •« 


Wh»?n  m  constraints  are  applied  on  the  ,n,  dependent  parameters, 

X,  we  reduce  the  number  of  variables  to  k  (n-m)  independent  parameters 
V.  The  normal  equations  are  now 

A,'.v  ••  B 


A  and  B  must  be  recalculated  in  terms  of  the  original  derivative 
matrix  D.  The  new  derivative  matrix  E,  in  terms  of  the  K  independent 
variable  v  is  given  by 


where  is  a  transformation  matrix  relating  the  derivatives  of  the 
dependent  parameters  with  respect  to  the  independent  parameters. 

The  derivative  of  E  with  respect  to  the  independent  parameters  is 
obtained  by  the  chain  rule. 


A 


eS: 


B  -  iF  E*^ 


The  dependent  solutions  are  given  by 


Ax  »  JAv 


The  parameter  change  vector  can  be  obtained  by  matrix  inv’ersion 
or  the  eigenvector-eigenvalue  method 

1.  Matrix  Inversion  (INVERT) 

A  Ax  "  B 

(NAN)  (n"^  ax)  -  N  •  B 


AX 


A  X 


(N/.N) 


AX  •  N  •  ax' 


S6 


•  (NB) 


2.  Eiyenvector-eiij'^f'Vdlue  (IGN)  : 


A  Ax  ■  B 


-1 


(NAN)  (N  Ax)  «•  NB 


T"^  (NAN)  T  (n"^*-X)  -  T 


\  -  (NAN)  T 


^-1 

T  -  T 


V  •  (n"^  Ax)  -  { 


Ax"  -  (n"^  Ax)  -  z  •  x"^ 


Ax  ■  N  T  x" 


M 

o  -  (  AF^/M-n  ♦  ml 
1-1 


Standard  deviation  of  the  parameter 


S 


j 


r 


Correlation  matrix 


VIII. 


Proqram  packaqe  consists  of: 

1.  This  users  manual. 

2.  A  magnetic  tape  with  Fortran  codinq  and  input/output 
files  for  an  electron  diffraction  analysis  of  silicon 
and  an  x-ray  diffraction  analysis  of  silica  glass. 

2.  A  printed  listing  of  these  files. 
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